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An Experimental Study of Atmospheric Homogeneous Nucleation; Cluster 
Growth and Gas-Particle Reactions of H2SO4 

Final Technical Report 

The work proposed on this project included both field and laboratory studies. The 
laboratory studies were to consist of measurements of H2SO4 uptake and evaporation from 

aerosols of varying chemical composition, while the field component would include 
measurements of H2SO4 and other compounds which would be conducted as part of a large field 

campaign. The exact field campaign in which we would participate could not be identified at the 
time the proposal was written, but it was hoped that participation in a larger campaign would 
provide additional knowledge of the chemical and physical environment in which the 
^SCVaerosol interactions under study were taking place. All studies proposed would also be 

conducted jointly with Dr. Peter McMurry’s aerosol group. By chance, the opportunity to 
conduct such an H 2 S0 4 /aerosol/ultrafine particle study in conjunction with an OH 

intercomparison/photochemistry study became available very early in this project (September 
1993 ). This study was conducted at Caribou, Colorado in conjunction with several other groups 
from NCAR, NOAA and a number of universities. Our group measured OH, H2SO4, SO2, and 
H2O, while Dr. McMurry’s group measured ultrafine particles, and total particle number and size 
distribution. In addition measurements of HO2/RO2, 03, NO, NO2, NO y , CO, hydrocarbons, 
CH 2 0 , and other chemical compounds and meteorological parameters were performed by the 
other participants and a new laser oblation/mass spectrometry technique was also employed by 
the NOAA Aeronomy Laboratory to study aerosol composition. The study of aerosol production 
and growth in conjunction with photochemical measurements is highly advantageous because 
particle growth precursors such as H 2 S 0 4 or MSA are formed by OH initiated sulfur oxidation. 

The large number of hydrocarbon measurements included in this study were also important in 
understanding particle growth as discussed in Appendix A. 

Participation in this study produced a wealth of new understanding in a variety of aerosol 
growth related areas, each of which is summarized below: 

Balancing the production and loss of H2SO 

The major sources and sinks for gas phase H 2 S 0 4 in a remote measurement site are 
shown below: 
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step 
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H 2 S0 4 + ? ^ >artlC — — new particle 
nucleation 
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Reaction sequence one is reasonably well understood and the SO 2 + OH reaction is the rate 
limiting step. Reaction three is not well understood, but its net result is at least quantifiable in 
terms of an accommodation coefficient and a particle distribution. Reaction sequence two and 
four, however, are particularly poorly characterized. Thus, an attempt to balance all sources and 
sinks of H 2 S0 4 would prove most difficult. At a remote continental site, however, the 
concentration of S0 2 is typically much larger than that of DMS and therefore reaction sequence 
two is expected to make only a minor contribution. Continental air masses are also typically 
expected to have aerosol concentrations sufficiently high to make H 2 S0 4 loss through reaction 
sequence three much larger than through four. As will be shown in the following section, 
particle nucleation can occur in continental air masses, but the ultrafine particles formed have 
much lower concentrations and surface areas than the residual particles. Thus, in a remote 
continental site, the production of H 2 S0 4 through reaction sequence one should be 
approximately equal to its loss onto existing particles (3). The production of H 2 S0 4 can then be 
equated with the loss rate for steady state conditions: 

[0H][S0 2 ]kj = [H 2 S0 4 ]k3 
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where kj is the reaction rate coefficient of SO2 + OH and is the rate of loss onto particles. 

The ability to simultaneously measure all of the above parameters i.e., OH, SO2, H2SO4 and 
particle number and size distribution has made a test of our understanding of H2SO4 production 
and loss possible. Also, fortuitously OH and H2SO4 concentrations are found on the opposite 
side of the above equation. These two compounds are measured using the same instrument and 
reaction scheme, therefore, the relative uncertainties introduced into the above equality by these 
two measurements is quite small. 1 Unfortunately, however, the SO2 concentration and the 
aerosol number density and size distribution measurements probably still introduce about a factor 
of 2 uncertainty into this calculation, particularly at low SO2 concentrations (a commercial 
Thermal Electron 43 S Instrument with a detection limit of 100 pptv was used to measure S0 2 
and particle size and number density was measured with a University of Minnesota DMA). 
Wi thin a factor of 2, however, the following tests provided the first real indication of our 
understanding of H2SO4 production and loss. 

Figure la shows a comparison between measured H2SO4 concentrations and those 
calculated using measured OH and SO2 concentrations, aerosol size and number distribution, and 
which assume an accommodation coefficient of 1 .0. These same measurements were also 
performed at Mauna Loa Observatory in Hawaii in 1992; these results are plotted in figure lb in 
the same format used in figure la. In both cases, agreement is often seen to be better than a 
factor of 2. These measurements are described in more detail in Appendix B. 

The relatively good overall agreement both in the variability and absolute value of 
calculated H 2 S0 4 production and loss provides hope that this same approach might be applicable 
in estimating H 2 S0 4 production from DMS (if S0 2 sources can be accounted for using the above 
method). Data from the recently completed ACE-I campaign presently provide an ideal 
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opportunity to test these ideas, offering far more sensitive S0 2 and DMS measurements, multiple 
particle number and size distribution measurements, and more sensitive OH and H2SO4 
measurements. 

Particle Nucleation 

While the previous discussion of H2SO4 uptake by particles tends to support classical 
understanding of particle growth, recent nucleation studies are suggesting mechanisms quite 
different from classic bimolecular nucleation theory. By making simultaneous measurements of 
gas phase H2SO4 and ultrafine particles, it was quickly determined that, as expected, these two 
parameters were closely correlated. The concentration of ultrafine particles, however, depends 
both on production (as shown by reaction sequence four) and loss onto existing particle surfaces. 
Therefore it is not the concentration of ultrafine particles but rather only their rate of production 
that should depend on H 2 S0 4 concentration directly. 

Figures 2a and b show the calculated flux (j) of ultrafine particles passing through the 3 
to 4 nm diameter size range as a function of the H 2 S0 4 concentration at the Mauna Loa and 
Colorado measurement sites respectively. Both plots show a compilation of all relevant data 
obtained during the campaigns, and j was calculated from the measured particle distribution 
function and by assuming an H2SO4 accommodation coefficient of 1.0 for ultrafine growth. 

Since 3-4 nm is the smallest size range of particles that can be measured, the flux through this 
size range presently provides the best estimate of nucleation rate. Also, however, because 
nucleating particles can be lost onto varying concentrations of larger particles prior to being 
detected in the 3-4 nm range, it is assumed that the highest calculated fluxes best characterize the 
actual nucleation rate (see Appendix C). If a line is drawn approximately through these high flux 
points in both figures 2a and b, its slope in both cases appears to be about 2. This suggests a 
nucleation rate which varies approximately as the square of the H2SO4 concentration. This is in 
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sharp contrast to the far more vertical lines showing classical nucleation theory predictions. 
Present results suggest a far weaker functional dependence on H 2 S 0 4 and also that nucleation 
can occur at much lower absolute concentration of H2SO4 than predicted by classical theory. If 
an accommodation coefficient of 0.5 were assumed in the calculation of j, the slope of the 
maximum flux lines would remain the same (2) and the j values would only drop by a factor of 
two. These calculations are discussed in more detail in Appendix C. 

The large disparity between classical bimolecular nucleation theory and the present 
results requires a re-evaluation of the use of this theory for describing nucleation in the 
troposphere or at least the lower troposphere where the present measurements were performed. 
To pose this issue in a different way: why, in hindsight would molecular clusters of H2SO4 and 
H 2 0 not be expected to react with basic compounds such as NH3? NH3 is typically more 
abundant than gas phase H 2 S 0 4 , it is highly soluble, and is commonly found along with H 2 S 0 4 
in larger particles. If NH 3 were to react with a molecular cluster containing nucleation 
precursors such as H2SO4 and H 2 0 , it could form ammonium bisulfate or even ammonium 
sulfate. This might well reduce the vapor pressure or the loss rate of H2SO4 from these clusters. 
If this were to occur, the ammonia could have a significant effect on the nucleation process such 
as that shown in figure 2 . This possibility is discussed in more detail in the laboratory results 
section. 

Ultrafine Growth Rates 

Small particles (particle diameter « mean free path) are typically assumed to grow at a 
constant diameter growth rate, independent of their size; i.e., they incorporate molecules at a rate 
which is proportional to their surface area. The most conspicuous growth thus occurs when 
aerosols are the smallest. If nucleation is assumed to begin with a cluster of just a few H 2 S 0 4 , 
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H2O, and maybe some NH3 molecules, then the initial size of this nucleating particle is ~ 1 nm. 

It need then only grow an additional 2 nm in diameter before it can be detected as a particle. 
Small ultrafine particles therefore offer a unique opportunity to study early stages of growth. 
Since the nucleation process requires the presence of relatively high H2SO4 concentrations, it is 
expected to have at least as strong a diurnal variation as H2SO4 and probably stronger 
([H2SO4] 2 ). Thus, the time between sunrise and the appearance of the first ultrafines provides a 
measure of the approximate time required for growth into the 3-4 nm range or ~ 2 nm of growth, 
under the ambient conditions for which the observations were made. Figure 3 a shows a 
compilation of all the H2SO4 data and the average clear sky solar fluxes plotted as a function of 
time of day for the previously discussed 1993 Colorado measurements. Note the close 
correspondence between the rise of solar flux and H2SO4 concentration. Unfortunately, 
afternoon and evening were often cloudy and/or rainy and therefore H2SO4 values did not 
decrease with clear sky solar flux. Figure 3 b shows a similar plot for ultrafine particles. These 
do not begin to rise at sunrise with H 2 S 0 4 . Rather, they change little until one or two hours after 
H2SO4 has risen. Using this rise time, a 1-2 nm diameter/hour growth rate is estimated for 
ultrafine particles at this site. Maybe a slightly better estimate of growth rate would be obtained 
by assuming that on average the nucleation process begins when the average H2SO4 
concentrations shown in figure 2 a were reached. This would suggest growth rates closer to 2 
nm/hour. Interestingly, this rapid growth rate can not be obtained by the incorporation of 
H 2 S 0 4 , H 2 0 or even NH 3 alone. This observed growth rate is about an order of magnitude 
larger than can be explained by H2SO4 incorporation, even with an accommodation coefficient 
of one. Though the data were more scattered at the MLOPEX II study and probably more 
dominated by transport, this large disparity between observed and calculated H2SO4 induced 
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growth did not arise. At present, the most likely explanation for this very rapid observed growth 
rate is that other low volatility compounds are contributing to the growth process. Since this 
disparity appears very large only in continental air masses, very large or highly polar 
hydrocarbons would seem likely candidates to account for this rapid growth. This possibility is 
discussed in more detail in Appendix A. 

In November and December of 1995 the ACE-I (Aerosol Characterization Experiment) 
was conducted in a remote southern hemisphere marine environment with many goals in 
common with the present project. After our past success it was decided that participation in this 
study could be even more productive than the past study particularly since it also provide an 
opportunity to contrast continental vs marine environments. Our instrument made measurement 
of OH, H 2 SO 4 , and MSA for nearly 300 hours on the NCAR C- 130 as part of this field 
campaign. While our measurements were highly successful much of the data from other 
investigators is not yet available. Therefore, most of the data interpretation is not yet possible. 
Figure 4a, b, and c, however, show some typical measurement results for OH, H 2 SO 4 , and MSA 

respectively. 

Laboratory studies in support of and sometime even helping to guide field measurement 
strategies (such as the inclusion of NH 3 measurements in ACE-I) have also been highly 
successful. One of the major questions raised by the 1993 campaign and explored by ACE-I, is 
the possible role of NH 3 in particle nucleation. To investigate this possibility, joint laboratory 
studies were conducted to measure the vapor pressure of H 2 SO 4 above aerosols containing 
various mixtures of H 2 S0 4 and NH 3 . Figure 5 shows the results of these studies for several 
different ratios of NH 3 to H 2 S0 4 from 0-2. In each case, the H 2 SO 4 vapor pressure is plotted as 
a function of relative humidity. For pure H 2 SO 4 the vapor pressure is compared to that 
calculated by Ayer et al ., 2 and agreement is seen to be reasonably good. As the ammonium 
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content of the particles is increased to 13, 15 or 33%, the vapor pressure of H 2 SO 4 appears to 
drop slightly. However, when ammonia concentration reaches 80% of H 2 SO 4 , the drop becomes 
very pronounced. When the average composition of ammonium bisulfate is reached, the vapor 
pressure of H 2 SO 4 is reduced by two-three orders of magnitude, depending on relative humidity. 
To a crude approximation, this dependence on ammonia content suggests that each ammonia 
molecule dramatically reduces the ability of one sulfuric acid molecule to leave the particle with 
the excess acid molecules loss largely unimpeded. At still higher ammonium/sulfate ratios (2) 
the H 2 SO 4 vapor pressure drops so low that it is difficult to measure. While the H 2 SO 4 vapor 
pressure above ammonium sulfate is very low, the vapor pressure of ammonia is probably quite 
high . 3 It is therefore less likely that ambient NH 3 concentrations will be sufficiently high such 
that ammonium sulfate molecular clusters will play an important role in atmospheric nucleation 
(though little quantitative can be said about ammonia vapor pressures until additional 
measurements are made). The two to three order of magnitude reduction in H 2 SO 4 vapor 
pressures for the ratio of 1 : 1 shown in figure 5 is, however, already sufficient to dramatically 
alter calculated nucleation rates if this same vapor pressure dependence applies to molecular 
clusters/ultrafine particles. 

In attempting to understand how NH 3 stabilized molecular clusters might grow into 
particles, one must also briefly examine the role that statistics might play in this process. For the 
large particles shown in figure 3 (10 2 nm diameter) ammonium content is listed in fractional 
concentrations. However, for small molecular clusters with only a few H 2 SO 4 molecules, there 
is presumably a probability distribution for having 1 , 2 , 3, — ammonia molecules associated 
with each cluster. If the molecular cluster of interest has no ammonia molecules or perhaps only 
one ammonia molecule compared to 2 or 3 sulfuric acid molecules, then from the predictions of 
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classical theory shown in figure 2, it probably has essentially no possibility of growth for the 
H2SO4 concentrations typically measured. However, even if the average NH3/H2SO4 
concentration ratio is well below one there will still be some significant fraction of molecular 
clusters containing, for example, two H2SO4 molecules and two NH3 molecules. If the results 
shown in figure 5 can be applied to molecular cluster growth, a reduction in H2SO4 loss by 
several orders of magnitude might well make the further growth of the latter cluster probable. 
While highly speculative, a scenario such as this could explain the present results, which suggest 
a kinetically controlled process depending on [H2SO4] 2 and, at the same time, predict a 
nucleation rate which is many orders of magnitude slower than that calculated for H2SO4, uptake 
by the clusters at the collision rate. It is also interesting to note that the calculated j values in 
figure 2b are somewhat higher than in figure 2a, and the ammonia concentrations are also 
expected to be much higher over the continent than at Mauna Loa. 4 ’ 5 These measurements are 
discussed in more detail in Appendix D. 

Most recently, measurements of the uptake of H2SO4 on sodium chloride and ammonium 
sulfate particles have been conducted. The preliminary results of these studies suggest an 
accommodation coefficient in the 0 . 5 - 1 .0 range which is in good agreement with the previous 
field study results for balancing the production and loss of H2SO4 (see Appendix B). Analysis 
of these data are still underway but the results will be written up and published once the analysis 
is completed. 

Conclusions 

Recently developed instrumentation is opening up exciting new avenues for investigating 
particle nucleation and growth. In the absence of significant DMS and particle nucleation, 
H2SO4 production and loss can be reasonably well-characterized and balanced. This offers 
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hope that the production rate of H 2 S0 4 from DMS might also be quantifiable in the near future. 
The application of classical nucleation theory in the troposphere has been questioned, and a more 
kinp.tir.ally controlled process possibly involving NH3 is posed as a possible alternative. The 
results of ACE-I which included NH3 measurements along with those of H2SO4, H 2 0 and 
ultrafines are expected to substantially improve present understanding in this area. If the results 
shown in figure 2 represent tropospheric nucleation in general, then nucleation events may be far 
more prevalent and important than present models would suggest. Finally, the measurement of 

* 

small, size selected, ultrafine particles now opens up the possibility of studying particle growth 
processes on a time scale which is short compared to a day. This is extremely important, 
because particle growth precursors and dynamics both vary substantially over the time period of 
a day. Thus, attempting to understand an already highly complex process with multiple varying 
parameters is particularly difficult. The ability to determine growth over the period of an hour or 
two under relatively stable chemical and dynamic conditions, however, offers far more hope of 
obtaining a detailed understanding of this process over the next decade. 

A great deal of new insight into particle nucleation and growth has been gained by this 
group and others over the past half decade and far more progress is anticipated in the near future. 
Some of the areas in which this progress is desired and anticipated are: understanding DMS 
oxidation and the subsequent yield of H 2 S0 4 and MSA; quantifying the role of NH 3 in particle 
nucleation; being able to accurately predict nucleation events from precursors or even source gas 
measurements; and understanding and quantifying the role of hydrocarbons in aerosol growth. 
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Abstract. Ultrafine aerosols, with diameters below 10 nm, nucleate from gas phase species. The 
composition of newly formed ultrafine atmospheric aerosols is not known with certainty; new 
particles have variously been conjectured to be sulfates, organic compounds, and sulfate/organic 
mixtures. The 1993 Tropospheric OH Photochemistry Experiment at Idaho Hill, Colorado 
provided an opportunity to examine the question of which class of compounds— sulfates or 
organics— make the major contribution to new particle formation in the unpolluted troposphere. 
This study compared the production rates of sulfuric acid (from the oxidation of sulfur dioxide) 
and oxidized organic compounds to gauge their relative contributions to the formation of ultrafine 
particles. Potential organic precursor species examined in this study were the naturally occurring 
terpenes a- and p pinene, and the anthropogenic hydrocarbons toluene, m-xylene, ethyl benzene, 
1,2,4 trimethyl benzene and methylcyclohexane. The calculated production of oxidized organics 
appeared well correlated with total particle surface area and volume, suggesting that at least some 
of the organic compounds formed in gas phase reactions condensed upon the preexisting aerosol. 
New particle formation was found to be more highly associated with elevated production of gas 
phase sulfuric acidfvia the SO 2 -OH reaction) than with production of oxidized organic products, 

although data from one day, during which sulfuric acid production and total aerosol surface area 
were both lower than usual, provided evidence for the involvement of terpene species in new 
particle formation. The results suggest that for this continental site, sulfuric acid was probably 
responsible for most of the observed new ultrafine particle formation. Low volatility organic 
compounds may have caused particle formation under the right conditions, but were more likely to 
condense upon pre-existing particles. 



Introduction 


Atmospheric aerosols play a significant role in shaping the Earth's climate, both by the 
direct scattering of solar energy [Charlson et al., 1992] and through their role as cloud 
condensation nuclei [Twomey et al., 1984]. Since atmospheric particles are continually subject to 
coagulation, scavenging by larger particles, and removal by precipitation, a source of new aerosols 
is evidently required to maintain observed particle size distributions in the background troposphere. 
These newly formed particles comprise the ultrafine aerosol. With diameters less than 10 nm, 
ultrafine particles lie at the small size extreme of the atmospheric aerosol. They are formed from 
gas phase precursors, having grown past the critical cluster size to attain stability as particles. 

Since ultrafine particles are quickly lost by diffusion to the surfaces of larger particles, their 
presence in regions far from primary emissions is evidence that nucleation from gas phase 
precursor species is taking place. 

The sources and composition of the ultrafine aerosol have been subject to much debate. 

The question attracted attention early in the history of aerosol science. In the mid nineteenth 
century, John Tyndall passed intense light through a vessel of organic vapors and observed the 
formation of a light blue cloud, which gradually turned whitish over time [Tyndall, 1868]. He was 
observing the formation of small particles which scattered strongly in the short wavelengths, 
followed by particle growth to the Mie scattering regime. This arcane experiment became known 
to milli ons living in the Los Angeles basin when Haagen-Smit and Fox [1956] established the 
mechanisms for the various smog reactions, driven by the photochemistry of oxides of nitrogen 
and reactive hydrocarbons. The production of small particles in mixtures of organic vapors and 
oxidizing agents has been repeated in the laboratory many times in the last four decades. Went 
[1960] proposed that the same reactions may produce sub-micrometer aerosols by oxidation of 
terpenes, naturally occurring hydrocarbons emitted by many species of vegetation. He and others 
[Rasmussen and Went, 1965; Went, 1966; Lopez et al., 1985] have argued that these reactions are 
the cause of the small particles comprising the "blue haze", common over forested areas. 



Laboratory studies of hydrocarbon reactions with the hydroxyl radical (OH) and/or ozone 
(O3) have observed the production of condensable species and aerosols. Naturally occurring 
monoterpenes (C10H16) and sesquiterpenes (C15H24) may be oxidized by OH or O3 to form a 
large number of terpenoid products [Graedel, 1979; Yokouchi and Ambe, 1985; Hatakeyama et 
al., 1991; Pandis et al., 1991; Palen et al., 1992; Grosjean et al., 1992; Shu and Atkinson, 1995]. 
In addition, OH and ozone reactions with many anthropogenic hydrocarbons have been observed 
to form aerosols [Stem et al, 1987; Wang et al, 1992]. The products of these organic reactions 
possess a wide range of vapor pressures, making it difficult to specify the condensing compound. 
Two probable reaction products of a-pinene and OH , pinonic acid and pinonaldehyde, have been 
identified in atmospheric aerosols [Wilson et al., 1972; Hatakeyama et al., 1989]. 

In the 1970s attention within the aerosol community generally shifted from terpenes to the 

sulfate aerosol. New ultrafine particles in the background atmosphere were postulated to nucleate 
from vapor phase water and sulfuric acid (H^SO^), the latter formed by the oxidation of sulfur 

dioxide (S0 2 ) by OH [Doyle, 1961; Mirabel and Katz, 1974; Jaecker-Voirol, 1989], Subsequent 

reaction with gas phase ammonia results pi a partially to totally neutralized particle, the latter 
usually as ammonium sulfate, (NH^SO^. Most recent modeling studies of new particle 

formation in the atmosphere have used modifications of classical binary nucleation theory 
involving water and sulfuric acid [Kreidenweis et al., 1991; Raes and Van Dingenen, 1992; Easter 
and Peters, 1994; Raes, 1995]. 

No direct composition information exists for ultrafine particles. Composition 
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measurements of particles down to the lower practical size limit (about 0.05 Jim to date) suggest 
that submicron particles have mixed compositions, with orgaruc/sulfate ratios that may vary 
according to particle size and sampling location. Weiss and co-workers [1977] examined aerosols 
from three rural locations in the southern and midwestem United States. Their method focused on 
particles between 0.1 and 1.0 pm in diameter, which includes the bulk of the optically important 
haze producing aerosol. The authors found that sulfate dominated the sub-micrometer aerosol 
under all wind and meteorological conditions, consistent with an earlier study [Charlson et al., 



1974]. On the other hand, size resolved particle composition measurements in the western United 
States (the Los Angeles area and near the Grand Canyon) by Zhang and coworkers [1993] found 
that carbon comprises a significant fraction of the aerosol mass in the diameter range 0.05 to 1.0 
Jim. Earlier studies [Hoffman and Duce, 1974; Ketseridis and Eichmann, 1978] also showed 
considerable organic material in particles smaller than 1 pm. Particle composition was measured in 
a number of U.S. national parks and wilderness areas by Malm et al [1994]; the authors found the 
largest single component of aerosols under 2.5 pm in diameter to be sulfates in the eastern United 
States and organic compounds in the west. Recent measurements [Novakov and Penner, 1993; 
Rivera-Carpio et al, 1994] suggest that organic carbon may be a major constituent of cloud 
condensation nuclei. 

While composition data on particles large enough to be measured may yield insight about 

particle growth, they tell us little about the nucleation phenomena that form new particles from the 

gas phase. In the absence of such data for nanometer-sized particles, inferences about the 

composition of newly formed particles may be made through simultaneous measurements of 

ultrafine particles and possible gas phase precursors. The opportunity for such a study arose 

during the 1993 Tropospheric OH Photochemistry Experiment at Idaho Hill, Colorado. For most 

of the month of September 1993, concurrent measurements were made of ultrafine particles, fine 
particles (with diameters between 15 and 500 nm), gas phase S0 2 and I^SC^, ozone, OH, and a 

number of nonmethane hydrocarbons. Wind, temperature and other meteorological data were also 

collected. This paper analyzes the Idaho Hill data set to gain insights into the contributions of 

possible precursor gas reactions to new particle formation at this relatively clean continental site. 

Correlations between particle formation rates and steady state sulfuric acid concentrations at 

Idaho Hill were the subject of an extensive analysis by Weber et al [1996]. The study found that 

new particle formation was well correlated with the concentration of vapor phase sulfuric acid. 
However, the study also found 1) particle nucleation occurred at I^SO^ concentrations much 

lower than predicted by classical binary nucleation theory for the H^SO^ -H 2 O system; 2) particle 

nucleation rate appeared to depend on the square of H^SO^ concentrations, rather than the higher 



power indicated by theory; and 3) particle growth after nucleation was eight to thirteen times faster 
than could be explained by condensation of and its associated water. The authors 

concluded that new particles at Idaho Hill formed from sulfates, possibly with the involvement of 
other species, such as atmospheric ammonia. It was hypothesized that condensation of organics 
might contribute to the observed high growth rates of ultrafine particles, although organic species 
were not explicitly considered in the analysis. The present work extends the analysis of Weber et 
al, examining the contributions of sulfate and organic precursor species to new particle formation at 
Idaho Hill. 


Apparatus and Techniques 

The site . The measurements were made during September 1993 on a forested ridge in the 
Rocky Mountains about 25 km west of Boulder, Colorado at a pressure elevation of about 700 mb 
(3070m above mean sea level). Detailed descriptions of the physical and meteorological conditions 
at the site are found in Mount and Williams [this issue] and Olson et al [this issue]. 

The local topography tended to force ground level winds into either an easterly (upslope) or 
a westerly (downslope) direction. Areas to the west of the site were sparsely populated, so the 
westerly downslope winds brought clean, nearly background continental air to the site. The 
easterly upslope winds were associated with anthropogenic influences from the Denver metro area 
and from agricultural land within Boulder County. The arrival of upslope winds shows clearly in 
the aerosol record as an increase in the concentration of anthropogenic hydrocarbons and particles 
with diameters greater than 50 nm. 

Aerosol apparatus . The concentration of ultrafine (hereafter UF) aerosols was measured 
with an ultrafine condensation nucleus counter (UFCNC), a modified commercial CNC (Model 
3020, TSI, St. Paul, MN) which had a lower detection limit of about 2.7 nm (50% detection 
threshold of 3 nm) [Stolzenburg and McMurry, 1991]. Samples were collected with a 2.5 cm 
inlet about 5 m above ground level. The aerosol sample flow was taken from the centerline of the 
inlet flow to minimize the effects of diffusional losses to tubing walls. Measurements were made 



continuously with time resolution of about 2.5 minutes. The UFCNC was capable of detecting 
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UF particles (3 to 4 nm diameter) at a concentration as low as about 0. 1 cm . 

Like all CNCs, the UFCNC gives a count of total particles greater than a threshold size by 
condensing a working fluid (butyl alcohol in this case) on these particles; the resulting droplets are 
sensed optically. While CNCs allow optical detection of particles much smaller than the 
wavelength of light, they do so at the expense of size information on the original aerosol. The 
UFCNC partially overcomes this limitation of conventional CNCs by modifying the aerosol flow 
through the instrument's condenser region, confining all particles to the condenser centerline. This 
ensures that all particles are exposed to about the same peak saturation ratio and have similar 
transport times through the condenser. Along the centerline, input particles encounter a rising 
supersaturation that peaks about 2/3 of the length of the condenser. Smaller input particles, which 
require higher supersaturations to nucleate droplets, nucleate later than larger particles. This results 
in a final droplet size that reflects the size of the original input particle, allowing a method of UF 
aerosol size discrimination [Brockmann, 1981; Stolzenburg, 1988; Ahn and Liu, 1990]. In 
practice, particles with diameters larger than about 12 nm all grow to about the same final droplet 
size, but the sizes of particles smaller than this may be measured with the UFCNC by analyzing the 
height of the pulse produced at the optical detector [Saros et al, 1996; Weber et al., 1995a]. The 
optical output from the UFCNC was coupled to a multichannel analyzer (EG&G model 916A) for 
pulse height analysis. 

The response of the UFCNC to particles of sub- 15 nm diameter was determined using 
monodisperse test aerosols. The concentration of newly formed aerosol was calculated from the 
count rate in the portion of the pulse height spectrum corresponding to particles between about 3 
and 4 nm diameter. While this size is much larger than that of the critical molecular cluster which 
would truly represent a new nucleated particle, it is the nearest to the newly formed particle size 
range we can observe with our instrument. Throughout this study, measured concentrations of 3 
to 4 nm diameter particles are taken to be the concentration of newly formed aerosol (the time delay 
between actual nucleation and attainment of detectable size is addressed in Results and Discussion, 



below). Experimental uncertainty increased with decreasing particle concentration; the relative 
uncertainty for a concentration of 0.1 cm'3 was about 15% . 

Measurements of fine particle size distributions were made with a scanning electrical 
mobility spectrometer (SEMS), consisting of a differential mobility analyzer (DMA) and a CNC 
(TSI 3760). In the SEMS, the DMA is continuously scanned through its voltage range as the CNC 
detects particle concentrations in each of several discrete size ranges [Wang and Flagan, 1990]. A 
full scan through the DMA voltage range (up and down) took about 2.5 minutes. For most of the 
experiment, the SEMS measured particles in the diameter range between 0.015 and 0.5 pm; on 
occasion, high particle loading required altered instrument flows which set the upper limit of the 
detected size range at 0.26 pm. SEMS data were inverted by the method of Hagan and Alofs 
[1983] to yield the fine particle size distributions. This inversion technique accounts for both 
multiple charging and aerosol transport losses in the DMA. From these distributions were 
calculated the total aerosol surface area and total volume. Both these quantities were dominated by 
particles above 0.1 pm in diameter, even though smaller particles (0.05 to 0.1 pm) were much 
more numerous throughout the experiment Measurements of fine particles, including those over 1 
pm in size, are treated in more detail elsewhere [Murphy, this issue]. 

Gas Phase Sulfuric Acid. Sulfur Dioxide. OH and Ozone. The gas phase concentrations of 
H^SO^ and OH were measured every five minutes with a chemical ionization mass spectrometer 

5 -3 

(CIMS). Detection thresholds for both species were about 10 molecules cm . More details on 
the operation of the CIMS can be found in Tanner et al [this issue] and [Eisele, 1993]. 
Uncertainties in gas phase concentration measurements were a maximum of +/- 35%. Additional 
OH concentration data was available from the long path measurements [Mount et al, this issue]. 
Sulfur dioxide was measured with a commercial pulsed fluorescence instrument (Thermo 
Environmental Instruments Model 43S) at about ten minute intervals. Uncertainties were usually 
under 20%, but at times ranged over 100% due to problems with instrument zero drift Ozone 
concentrations were measured both by long path differential absorption spectroscopy and a 



commercial instrument. The results of the two methods agreed to within about 2% [Harder et al, 
this issue]. 

Nonmethane Hydrocarbons (NMHCs\ The concentrations of a number of NMHCs were 
measured hourly by gas chromatography with cryogenic trapping; details of the hydrocarbon 
sampling are found in [Goldan and Kuster, this issue]. Of particular interest among the measured 
organic species are the monoterpenes, a— and (3 pinene. These two compounds are often invoked 
as candidates for precursors of particle nucleation (other potential terpene precursors, such as d- 
limonene and the sesquiterpenes, were not measured at Idaho Hill). Also examined were the 
aromatics trimethyl benzene, ethyl benzene, toluene, and m-xylene and the alkane 
methylcyclohexane. These (presumably) anthropogenic NMHCs have all been observed to react 
with OH or other oxidants to yield measurable aerosol mass. During the early part of the 
experiment, data were collected only during the hours 10 a.m. to 3 p.m.; sampling was expanded 
to 24 hour during the latter half of the field campaign. 

Analysis. Theoretical and laboratory studies of particle nucleation in chemically reacting 
systems [McMurry and Friedlander, 1979; McMurry, 1980; McMurry, 1983; Kreidenweis et al., 

1991; Wang et al., 1992] have shown that the maximum concentration of new particles formed in 
such systems (N max ) depends on 1) the vapor pressure of the condensable compound(s) formed 

by gas phase reactions, 2) the rate of formation of these compounds, and 3) the total surface area 
of any preexisting aerosol. In the present work, we examine the Idaho Hill data set to determine 
whether or not the daily N ma x of UF particles exhibit dependence on the formation rate of the 

nucleating species and total aerosol surface area. The aim of this investigation is to compare the 
relative importance to new particle formation of gas phase S02-OH reactions, which ultimately 

form sulfate particles, versus the oxidation reactions involving natural and anthropogenic NMHCs, 
a subset of which may form new particles. 

The peak value of UF particle concentration Nmax was determined for each of the 17 days 
for which complete data existed. Reaction rates R were calculated separately for each potential 
precursor-oxidant combination 



^PRECURSOR ~ ^ [precursor] [oxidant] 


using literature values for the rate coefficients k. Thus the Rso2-OH is given as 

R S02-OH = k l [S02][OH] (1) 

where kj is that given by [Seinfeld, 1986]. Sulfur dioxide-ozone reactions were not considered in 
this study. 

Unlike the sulfate system, reactions of the NMHCs do not yield a single condensable 
product To facilitate comparison between reacting systems, reaction rates were summed for the 
terpenes and the anthropogenic hydrocarbons. An aggregate formation rate for terpenoid products, 
Rterp* is defined as 

Rterp = Zn t0i kJ ti (2) 

'•* \ - 

where n t are the gas phase concentrations of the monoterpenes (a- or P-pinene), Oi the 
concentrations of OH or ozone, and kti the appropriate rate constant for each reaction. The terms 
fti are fractional aerosol coefficients (FACs), introduced as weighting factors to account for the 
varying tendency of these reactions to form aerosol products. Experimental values of FACs for the 
pinenes are those of Grosjean and Seinfeld [1989], Hatakeyama et al [1989], Arey et al [1990], 
and Grosjean et al [1992]. Similarly, an aggregate rate for the formation of condensable species 
from the anthropogenic hydrocarbons is defined as 

RaNTHRO = ^, n aQ ^ai/ai (3) 

aj 

where the na are concentrations of toluene, m-xylene, ethyl benzene, methylcyclohexane and 1,2,4 
trimethyl benzene. Rate constants and FACs are from Wang et al [1992], Stem et al [1987], and 



the above references. The rate constant and aerosol fraction for 1,2,4 trimethyl benzene were 
approximated using laboratory data for its 1,3,5 isomer. 

Mean values of the concentration of total aerosol surface area (in cm 2 cm -3 ) and volume (in 
cm 3 cm' 3 ) were also determined for each day of the experiment. When and if condensable species 
are produced in the gas phase, they may form new particles or condense upon preexisting particle 
surfaces. The latter is energetically favored, so new particle formation should be observed more 
often during periods of low total particle surface area, and the dependence of Nmax on R should be 
most pronounced under these conditions. Accordingly, the Nmax data for each reacting system 
(i.e., S0 2 -OH, terpenes, and anthropogenic hydrocarbons) were stratified according to the total 

surface area present during the measurement. Evidence of R dependence was taken to be an 
indication that the reacting species was involved in new particle formation. In this way, the relative 
roles played by organic and sulfate species in new particle formation may be directly compared. 

Results and Discussion 

A representative plot of UF partible concentration versus time of day is shown in Figure la. 
The Figure includes periods of downslope winds (Sept. 21, 24-25) and upslope wind conditions 
(Sept. 22-23). A diurnal variation is apparent in the UF data for the downslope days. Ultrafine 
particles typically appeared shortly after sunrise, increased to a peak between 8 and 10 a.m. local 
time, then decreased sharply after noon, usually vanishing by about 3 p.m. Such variation was 
observed on 80% of all experiment days; the portion rises to 85% if only downslope wind 
conditions are considered (9 of the total jS'days analyzed). Diurnal variations similar in shape, 
although delayed in time, have been observed with Aitkin nuclei in continental boundary layer air 
[Hogan, 1968; Marti, 1990]. The particle counter used in the present work is capable of detecting 
aerosols much smaller and at lower concentrations than those used in previous studies, which may 
explain the earlier daily onset and disappearance of measured UF particles. Ultrafine particle 
concentrations varied widely during upslope wind conditions, with the size distributions typically 
dominated by larger particles (diameter > 50 nm). 
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The data were examined for temporal relationships between UF concentrations and the R 
for each reacting system. Diurnal patterns in R s02 > ^terp * R^NTHRO would be expected 

due to the photochemical creation of OH and the diurnal changes in terpene emissions. If a given 
precursor is contributing to UF particle production, a time lag should be observed between 
increases in R and the first detection of UF particles, i.e., the time required for a newly nucleated 
particle to grow to the 3 nm lower detection limit. This growth time is a function of monomer size 
and gas phase concentration of the condensable species. For example, growth time for H 2 SO 4 

aerosols under conditions observed at Idaho Hill is estimated to be at least two to four hours, using 
the assumptions of collision controlled nucleation (for a discussion of collision controlled 
nucleation, the reader is referred to McMurry [1983] and Weber et al [1996]). In Figure la the 
daily variation of Rso2-OH IS coplotted with UF concentrations for the period of September 21-26. 
The SO 2 -OH reaction rates show a diurnal pattern. On seven of the 17 days analyzed, daily 
increases in R S q 2 -OH P receded those of UF particles, with a mean lead time of 1.9 hours. Figure 
lb shows daily variation in UF particles and R^rr Daily increases in Rjerp preceded those of 

UF particle concentration on five of 17 spiple days (by an average of 2.1 hours) and lagged UF 
concentration increases on 8 days. Figure lc compares the time variation of Ranthro U^ 

concentrations for the same period. Increases in the former preceded the latter on only one day 
during the experiment. On the basis of the data in figures la-c, terpenes appear to be the more 

likely than the anthropogenic NMHCs to be associated with new particle and growth. 

As noted in the previous section, we are looking for a positive dependence of N max on R 

for the possible condensing species. In Figure 2a, N may is plotted against the daily mean value of 
Rterr. Data collected under both upslope and downslope wind conditions are included in the 
Figure. The data have been segregated into two groups according to the total aerosol surface area, 
low (below about 5.0 xlO ' 7 cm 2 cm* 3 ) and high (greater than 6.0 xlO -7 cm 2 cm‘ 3 ). If terpene- 

OH/ozone reactions were associated with the production of new UF particles, a dependence of 

N on R should show up most clearly in the low surface area data. The data in Figure 2a do not 
max r 

seem to show such dependence; the data collected under low aerosol surface area conditions span a 



wide range of N m ax for a given R^erp- A linear fit of N max yielded near zero or negative 
dependence on R'yerp* 

Similar treatments of Ranthro ^d RS02-OH are presented in Figures 2b and 2c, 

respectively. As in Figure 2a, the data have been segregated into low and high values of total 
surface area. The data in Figure 2b exhibit little evident dependence of N ma x on RanthrO’ On 
the other hand, Figure 2c appears to show a slight N m ax dependence on Rs02-0H among the low 
area subset of the data. Although the data is sparse, a linear fit of N max to Rs02-OH f° r the low 
surface area data reveals a positive dependence, while the Rs02-OH dependence for the high 
surface area Nmax data is near zero. The scatter and scarcity of the SC>2 data do not permit 
rigorous statistical analysis; nevertheless this difference in R dependence distinguishes the SO 2 - 
OH data from that of the organic systems. 

Unlike with the organic systems, measurements of the SO 2 -OH system included data on 
its condensable product, H2S0 4 . To extend the examination of the role played by sulfate in new 
particle formation, the daily mean measured concentration of gas phase H 2 S0 4 at Idaho Hill is 
plotted against Nmax in Figure 3. Only (fata collected on days with low surface area conditions are 
included. A positive correlation of Nmax with [H 2 S0 4 ] is evident from the Figure. 

Figures 2 and 3 suggest that the new particle nucleation observed at Idaho Hill may have 
been more likely caused by SO 2 -OH reactions (that ultimately formed F^SO^) than by reactions 

involving the terpenes or the anthropogenic hydrocarbons. This is consistent with the analysis of 

Idaho Hill data by Weber et al [1996], where the authors found a correlation between measured 
concentrations of and UF particle concentration. 

If the organic compounds are reacting with ozone and/or OH to form a vapor species of 
low volatility, the latter may be condensing upon pre-existing particles, adding to the total aerosol 
surface area and volume. To check this, reaction rates R were compared with daily mean values of 
total particle surface area and volume for each day during which downslope (clean) wind 
conditions prevailed. All quantities were averaged from 6 a.m. to the peak daily value of R. 

Figure 4a shows both total surface area and total volume plotted against Rterp- Th e Figure 



shows that under clean air conditions, increased terpenoid production rates were associated with 
increases in particle surface area and volume. A similar plot of Ranthro is shown in Figure 4b. 
Note the low values of RANTHRO present in this downslope data (compare Figure 2b). Since the 
anthropogenic compounds included in Ranthro may be considered markers for upslope 
(polluted) air, this downslope data set considers only the low end of the Ranthro data range. 
Wi thin this range, Ranthro also appears correlated with particle surface area and volume. Figure 
4c shows that increasing Rso2-OH * s associated with increasing total aerosol volume, albeit 

weakly; little association with surface area can be seen in the Figure. Figures 4a-c suggest that as 
organic products were created, they may have preferentially condensed on pre-existing particles, 
adding to total surface area and volume, rather than nucleating new particles. This view is 
consistent with the lack of positive R dependence by Nmax found for the organic reactions in 

Figures 2a and b. It also fits with the time lag between observations of new particles and increases 
seen on about half of the sample days: as the products from pinene-OH/ozone reactions 

condensed on preexisting particles and increased the total particle surface area, new particle 
nucleation was increasingly inhibited; thus as R^np peaked, UF particle concentration was 

typically dropping. An alternate explanation for the observations shown in Figure 4 is that organic 
compounds may have been advected into the area along with larger particles, which dominate the 
aerosol area and volume distributions. Most of this effect, however, has been screened out of the 
data in Figure 4 by removing days with upslope wind conditions, which typically brought sharp 
increases in larger particles and anthropogenic hydrocarbons. 

While in general, oxidized organic products appeared not to play a primary role in new 
particle formation, there were possible exceptions. September 26 was a day of clean downslope 
wind conditions with low total particle surface area, about 2 x 10 -7 cmW 3 (the mean value for 
downslope wind conditions during the experiment was about 5 x 10* 7 cm 2 cm' 3 ). Sulfuric acid 
concentration was near the lowest seen during the experiment, averaging about 6 x 10 5 cm' 3 over 
the daylight hours. Figure 5 shows that the peak UF particle concentration was relatively low (70 
cm' 3 ) but non zero; R S0 2-0H remained below its mean value for d 0 ™ 11510 !* wind da y s > md its 




variations did not appear to precede or coincide with those of UF particle concentration. The 

values on this day were comparable to those of surrounding days. Unlike on most other 

days however, changes in on Sept. 26 led rather than lagged the appearance of UF 

particles, as shown in Figure 5. This correspondence, along with the lack of such correspondence 
with Rs02-OH the low concentration of H 2 SO 4 during this time, suggests that pinene 

reactions may have contributed to new particle formation seen on September 26. Increases in 

were observed to lead changes in UF particles on four other days (September 20, 21, 25, 

and 28). However, Rs02-OH H2SO4 concentrations were well correlated with UF particle 

concentration on these days as well, leaving open the question of which species may have played 

the major role in particle formation and growth on the days in question. On all other days with 
downslope winds, changes in Rs02-OH ^ H 2 SO 4 concentrations more closely corresponded to 

changes in new particle concentration than those of Rj^p • 


Conclusion 

If organic particle formation is an important source of new particles in the troposphere, a 
forested site in the western U.S. such as Idaho Hill may have been an ideal place to measure it 
Organic aerosols found in marine environments are thought to be largely of continental origin 
[Cachier et al, 1986]. The eastern U.S. has many anthropogenic sources of SO 2 which 

presumably cause sulfates to dominate new UF particle formation. Locations with significant 
vegetation and well removed from marine and anthropogenic sulfur emissions should offer the best 
setting to assess the role in new particle formation played by organic compounds in the undisturbed 
troposphere. The data from Idaho Hill do not provide strong evidence that this contribution is a 
major one. 

The sample size in the present study is admittedly small, and the data show a high degree of 
scatter, which limit the ability to make definitive statements about the ultimate source of new 
particles at Idaho Hill. However, while subject to caveats, the above analysis offers a fairly 
consistent picture of new particle formation. This study, and related work by Weber et al, suggest 



that in general sulfates probably dominated the observed new particle formation at Idaho Hill, 
while organic compounds appeared to be predominately involved with particle growth. Terpenoid 
compounds may have shown evidence of contributing to new particle formation on one day when 
sulfate particle production was suppressed. However, the condensable species formed when Gl- 
and P pinene were oxidized appeared more likely to condense upon preexisting particle surfaces, 
adding to particle volume and surface area. These findings suggest a mechanism for the mixed 
sulfate-organic particles observed by several researchers over the past two decades, and provide 
insight into a process not at present directly observable. 
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Figure Captions 

Figure 1. Representative sample of UF particle concentration, plotted with reaction rates R 
from sulfate and NMHC systems (see text for definition of R). The diurnal variation in particle 

concentration was observed on most days. 

a. UF particle concentration and sulfate reaction rate, Rso2-OH- 

b. UF particle concentration and terpene reaction rate, Rterp- 

c. UF particle concentration and reaction rate for anthropogenic hydrocarbons, 

Ranthro- 

Figure 2. The peak measured concentration of UF particles for each day (Nmax) of the 
Idaho Hill field study, plotted against daily mean values of R for each reacting system. The data 
are split into subsets according to the value of total aerosol surface area (see text). 

a. N m ax versus Rjerp- 

b. N m ax versus Ranthro- 
c- Nmax versus RS02-OIJ- 

Figure 3. The peak measured concentration of UF particles for each day (Nmax) plotted 
against daily mean values of gas phase H 2 SO 4 concentration. Only data collected during days of 
low particle surface area are included. 

Figure 4. Mean daily values of total particle surface area and total particle volume, plotted 
against the mean reaction rates of each reacting system. 

a. Area and volume versus terpene reaction rate, Rterp- 

b. Area and volume versus anthropogenic hydrocarbon reaction rate, 

Ranthro- 

c. Area and volume versus sulfate reaction rate, RS02-0H- 

Figure 5. Concentration of UF particles as a function of time on September 26. The 
calculated reaction rates Rterp ^d Rs02-OH ^ or day are also shown. 
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Abstract 

S im ultaneous measurements of aerosols and their expected gas phase precursors were 
made at Idaho Hill, Colorado, a remote continental site. This study used apparatus and 
tfchniq v*** similar to those employed in an earlier study at the Mauna Loa Observatory, Hawaii 
(Weber et aL, J. Atmos. ScL 52:2242, 1995 ). New particle formation, identified by the presence 
of ultrafine particles (nominally 3 to 4 nm diameter), was commonly observed in downslope 
(westerly) air and was correlated with high sulfuric acid (H2SO4) concentrations, low relative 
humidity and low aerosol particle surface area concentrations. The data point to H2SO4 as a 
principle nucleation precursor species whereas water vapor (H2O) apparently played a minor 
role. Particle production was observed at H2SO4 concentrations that are well below predicted 
values for binary nucleation of H2O and H2SO4, suggesting that another species was involved. 

Particle growth rates were estimated from the data with two independent approaches and 
in both cases were -5 to 10 times higher than can be explained by condensation of H2SO4 and its 
associated water. This suggests that species in addition to H2SO4 are contributing ultrafine 
particle growth. 

Finally, calculated steady-state H2SO4 concentrations were found to be in good 
agreement with measured values if the mass accommodation coefficient for H2SO4 on aerosol 
surfaces was assumed equal to -1. 



Introduction 

Tropospheric aerosol particles influence global climate by scattering radiation directly 
[Taylor & Penner, 1994] and indirectly by altering the scattering characteristics of clouds 
[Twomey et aL, 1984], Because both effects depend on particle size, processes which influence 
particle distributions can influence climate. The formation of new particles from gas phase 
precursor species plays an important role in regulating aerosol populations; in the remote 
atmosphere it is thought to be the primary source for new particles. The growth rates of newly 
formed particles are also critical since new particles must grow by orders of magnitude to 
influence the earth's radiation balance. 

A field study at the Mauna Loa Observatory, Hawaii, provided insight into particle 
formation in the remote marine troposphere [Weber et aL, 1995a], This effort has been extended 
by measurements made at Idaho Hill, Colorado, a remote continental site. The approach in both 
cases has been to examine atmospheric new particle formation by measuring the expected 
nucleation gas phase precursor species with simultaneous measurements of ultrafine (nominally 
3 to 4 nm diameter) and fine (-15 nm to -0.5 |im diameter) aerosols. Together, these two studies 
provide insight into atmospheric new particle formation and also permit comparison of particle 
production at a remote marine and continental site. 

Due in part to the extent of the earth's oceans, marine aerosols are thought to have a 
significant influence on global climate [Charlson et aL, 1987]. These aerosols have received 
special attention in both modeling and atmospheric measurements. Currently, classical binary 
nucleation theory involving sulfuric acid and water (H 2 SO 4 /H 2 O) [Jaecker-Voirol and Mirabel, 
1989] is for predicting atmospheric particle formation rates. It has not been demonstrated, 
however, that remote tropospheric particles are formed by this mechanism. Measurements by 
various researchers have provided evidence of an in-si tu source for marine aerosols. Clarke 
[1993] observed evidence of particle production in the Pacific free troposphere and modeling by 
Roes and coworkers [1992, 1995] suggests that the free troposphere is the major source for 
marine cloud condensation nuclei (CCN)- On Pacific Ocean cruises Covert et aL [1995] 



consistently observed evidence for particle formation in regions of subsiding air masses but saw 
no evidence of new particle formation in the marine boundary layer itself. Evidence of particle 
production has also been detected in the outflow regions of precipitating marine cumulus clouds 
[Perry and Hobbs, 1994] and in clear marine boundary layer air [Hoppel et aL, 1994a]. A 
common feature to the marine sites where researchers have detected new particle formation was 
low aerosol surface area concentrations (typically less than -5 |im 2 cm* 3 ). 

Many of the phenomena observed for marine aerosols may also occur in clean continental 
air, although less attention has been given to these regions. Both gas phase species and aerosol 
concentrations are generally higher in continental than in marine environments. In addition, a 
greater variety of species, both natural and anthropogenic, may interact with aerosols making 
continental regions more complex chemical systems. For example, vegetation derived 
hydrocarbons may also contribute to particle formation [Lopez et aL, 1985] and growth. Based 
on measurements at a mountain site in Puerto Rico, Novakov and Penner [1993] speculated that 
the majority of particles less than 70 nra diameter were organic. 

There is also evidence for new particle formation in clean continental air. A unique 
characteristic of these observations is a persistent diurnal variation in aerosol concentrations. 

This diurnal nature has been observed in a wide variety of continental settings and suggests a 
photochemical aerosol source [Koutsenogii and Jaenicke , 1994; Marti, 1990; Hogan 1968; 
Bradbury and Meuron, 1938]. As with marine environments, there is evidence of particle 
formation in the continental free troposphere [J Hofmann , 1993] and evidence of possible 
nucleation events in the vicinity of continental cumulus clouds [Radke and Hobbs , 1991]. 

This work reports on measurements made at Idaho Hill, Colorado, from 
5 September 1993 to 29 September 1993. Measurements representative of clean continental 
conditions are presented and compared to the earlier findings from the Mauna Loa Observatory. 



Site 


Idaho Hill is a remote continental site located on the eastern side of the Rocky Mountains 
(105*34*32" latitude, 39*58*56" longitude) slighdy east of the continental divide and 
approximately 25 km west of Boulder. The site is situated just below the timber line at an 
average pressure elevation of 0.7 atmosphere (3070 m elevation). The local topography consists 
of generally higher elevations to the west and lower elevations to the east The nearest populated 
areas (Boulder, Denver, etc.) are at lower elevations, east and south-east of the site. The western 
sector is sparsely populated and the prevailing westerly (downslope) winds were associated with 
relatively clean dry air. Upslope winds were typically from the SE, the direction of likely 
pollution sources. A more detailed description of the site is given by Mount and Williams 
[1995]. 

Instrumentation 

Pertinent gas phase species that were measured in this study include sulfur dioxide (S02), 
the hydroxyl radical (OH), sulfuric acid (H 2 SO 4 ), and water (H 2 O). In addition, alpha and beta 
pinene, which are also aerosol precursors, were measured. An analysis of the role of the organic 
species in new particle formation is reported elsewhere [Marti et aL, 1995]. Here the focus is on 
the role of H 2 SO 4 . 

Concentrations of molecular gas phase OH and H 2 SO 4 were measured to an absolute 
uncertain ty of ±35% using a unique chemical ionization mass spectrometer [Eisele and Tanner , 
1993; Tanner and Eisele , 1995]. The instrument's lower detection limit was approximately 
10 4 molecules per cm3. Sulfur dioxide was measured by pulsed fluorescence (Model 43S, 
Thermo Environmental Instruments Inc. Hopkington, Mass.). This instrument had a lower 
detection limit of -100 pptv (parts per trillion by volume). For the latter part of the study (24 
September -29 September) hourly S02 zeros were performed. This permitted correcting for a 
temperature-dependent systematic offset error observed with this instrument. The validity of OH 



me as urements with this instrument was recently confirmed in a methods intercoraparison study 
[Mount and Williams, 1995]. 

The ultrafine aerosol is defined here as particles in the narrow size range from 2.7 nm, the 
lower detection limit of the ultrafine condensation particle counter (UCPC), to nominally 4 nm 
diameter. Measurements of the ultrafine aerosol concentration were made with the prototype of 
the TSI 3025 UCPC [Stolzenburg and McMurry, 1991]. Ultrafine particle concentrations were 
determined by measuring the UCPC photo detector pulse heights [Weber et aL, 1995a] with a 
multichannel analyzer (MCA). This method of acquiring ultrafine particle size information is 
based on experimental work which showed that, for particles smaller than -15 nm, the final 
droplet size (pulse height), after growth in the UCPC condenser, is uniquely related to the initial 
particle size [Brockmann, 1981; Ahn and Liu, 1990; Stolzenburg 1988; Soros et aL, 1995]. 

For this study, we diluted the sample aerosol prior to measuring the ultrafine 
concentration. Nevertheless, on rare occasions, significant particle coincidence (more then one 
particle in the UCPC optical scattering volume at a time) resulted in inaccurate sizing of the 
aerosols. Based on laboratory experiments [Soros etaL, 1995], these episodes were 
identified by MCA dead times exceeding 15% and the data were ignored. The UCPC sample 
flow rate was 0.5 cm 3 s' 1 and the typical sampling period was 2.5 minutes. Considering only 
those uncertainties associated with Poisson counting statistics, these settings result in an 
uncertainty of ±35% for a measured ultrafine concentration of 0.1 cm‘3. Uncertainties decrease 
for concentrations higher than this. 

Bec ause ultrafine particles of a given size produce a range of pulse heights, our method of 
measuring the ultrafine particle concentration does not count all particles in the ultrafine size 
range (i.e., -3 to 4 nm diameter). This results in our underestimating the concentrations of these 
particles, typically by -50%. This shortcoming, however, is far outweighed by the 
measurement’s high sensitivity for detecting ultrafine particle concentrations. 

Aerosol size dis tributions covering diameters between 15 and —500 nm were measured 
with a scanning mobility particle spectrometer [Wang and Flagon, 1990]. A complete scan 



through this size range was typically completed every 2.5 minutes. These data were inverted 
using the method of Hagen and Alofs [1983]; aerosol surface area concentrations were obtained 
by integrating over these distributions, assuming spherical particles. The concentrations of larger 
particle's, -0.5 |im up to 9 (im diameter, were measured by other researchers at the site. These 
larger contributed at most 3% to the total aerosol surface area concentration when 
downslope conditions prevailed, although there were episodes when these larger particles 
(>0.5 pjn diameter) made significant contributions to the total aerosol surface area when the 
wind was from the polluted sector (upslope). 

Theory: Calculation of Particle Growth and Formation Rates 

The d at a are used to calculate particle growth and formation rates in a manner similar to 
that performed earlier with the Mauna Loa measurements [Weber et al., 1995a, 1995b]. 

If particle growth is limited by condensation of H 2 SO 4 vapor, the ultrafine particle 

growth rate (i.e., free molecular regime) is [Friedlander, 1977]: 

dDp _ 2avi(pi -Pd) (1) 

dt (2rankT) 1/i 

where a is the H 2 SO 4 mass accommodation coefficient, T the temperature, k is Boltzmann's 
constant, and m and v/, the mass and volume respectively of the condensing vapor. The 
measured partial pressure of H 2 SO 4 is pi, and the equilibrium partial pressure, pd- We assume 
an accommodation coefficient of one and that evaporation of H 2 SO 4 from the droplet is small 
compared to condensation (i.e., pd — 0.0). We also assume that growth is from condensation of 
only hydrated H 2 SO 4 molecules which are in equilibrium with the H 2 O vapor phase at 10°C. 
The degree of hydration is calculated using the equilibrium data of Gmitro and Vermeulen 
[1963]. These assumptions provide an upper limit if growth is only by condensation of hydrated 
H 2 SO 4 vapor. Actual growth rates, however, could be higher if additional species were also 
involved, or if significant numbers of H 2 SO 4 clusters also condense. 



The rate at which particles grow through a certain size due to condensing vapor 
can be determined from the particle size distribution function and particle growth rate. 
The formation rate of nominally 3 nm particles is. 


J(3nm) = 


dN 

dDn 


3.5nm 



\ 

'3.5nm 


( 2 ) 


where dN/dDp is the measured aerosol size distribution function at 3.5 nm and dDp/dt is 

the diam eter growth rate for 3.5 nm particles determined with Eq (1). We estimate the 
average size distribution function in the vicinity of 3.5 nm from our measurement of the 
ultra fine, aerosol concentration by: 


dN 

dD 


P /3.5nm 


_ measured concentration m 3 to 4 nm diameter range 

> n — 

4 nm - 3 nm 


( 3 ) 


Thus, our estimates of the 3 nm particle formation rate require simultaneous 
measurements of the H 2 SO 4 vapor and ultrafine aerosol concentrations. 


Results 

Aerosol Concentrations Versus Wind Direction 

At Idaho will measured ultrafine aerosol concentrations were found to depend strongly 
on wind direction. Figure 1 shows as a function of angular wind direction the concentration of 
particles greater than —3 nm (all counts from the UCPC) and the ultrafine aerosol concentration. 
In contrast to the total aerosol concentrations (Fig. la), which do not show a strong dependence 
on wind direction, high ultrafine concentrations were observed predominantly in downslope air 
(wind from western sector. Fig. lb). The narrow angular band where the highest ultrafine 
concentrations were recorded was the most common downslope wind direction. Total aerosol 
SM rfqr/> area was on average higher and more variable in upslope air. During the daytime 



between 600 and 1800 hours Mountain Standard Tune (MST), the period when significant 
photocheraically derived new particle formation is expected, the mean aerosol surface area 
concentration was 85±57 pm^ cm'^ for upslope air compared to 49±18 |Xm^ cm*3 for downslope 
air In this case the surface areas were calculated from particle size distributions spanning 15 nm 
to 9 nm diameter. The higher aerosol surface areas associated with upslope air were likely the 
result of anthropogenic sources to the east and are a possible cause for the low ultrafine aerosol 
concentrations when the wind was from this sector. Since the focus of this work is on particle 
formation and growth in the clean continental troposphere, only downslope data are analyzed 
further. 

Diurnal Variation of Species Concentrations at Idaho Hill 

Analysis of All Downslope Flow Data: In Fig. 2 the concentrations of OH, S02, 
H 2 SO 4 , and ultrafine aerosol for periods of downslope flow are plotted as functions of the local 
time of day. Also included with each plot is the clear sky ultraviolet solar intensity (wavelength 
range: 290 to 385 nm). This can be interpreted as the envelope of UV solar intensity since 
interference from clouds will result in attenuated values. 

The rapid increase and subsequent decrease in OH concentrations corresponding with 
sunrise (600 MST) and sunset (1800 MST) shown in Fig. 2a clearly indicate its photochemical 
source. Variations in OH concentrations are likely a result, in part, of variations in UV intensity. 
Similar observations of the diurnal variation of OH concentrations at the Mauna Loa Observatory 
are reported in Tanner and Eisele [1995]. Unlike OH, SO 2 had no diurnal pattern (Fig. 2b). 

Sulfuric acid is thought to be produced primarily by the reaction of OH and SO 2 (support 
of this is provided in the following section on predicted H 2 SO 4 concentrations). The diurnal 
variation of OH is responsible for the rapid rise in H 2 SO 4 levels observed at sunrise, Fig. 2 c. 
Note, there was practically no delay between the rise in OH and the earliest appearance of 
elevated H 2 SO 4 concentrations at sunrise. Scatter in H 2 SO 4 concentrations is probably from 
scatter in both the OH and SO 2 , and also the variation in the aerosol surface area concentrations. 



The ultrafine aerosol concentration also varied diumally. Fig. 2d. Unlike the H 2 SO 4 
concentration, which increased at sunrise (-615 MST, Fig. 2c), the concentration of ultrafine 
particles did not begin to increase significantly until approximately 730 MST, more than an hour 
la fpr We interpret this delay as the time required for the precursor vapor, H 2 SO 4 , (-1 nm 
diameter), to grow to the ultrafine size range (-3 nm diameter). In the following section this 
delay is use to estimate ultrafine particle growth rates. 

By la te afternoon the ultrafine concentrations began to decline and were low by nightfall. 
The pridian daytime ultrafine concentration was 48.2 cm'3 compared to 0.06 cm'3 during the 
night These highly mobile ultrafine particles have short life spans because of scavenging by 
existing aerosol surface. The median surface areas during day and night for downslope 
conditions were comparable; 42 and 48 jim^ cra'^ respectively. If there are no sources and the 
primary sink for the ultrafine aerosol is scavenging by preexisting particles, based on kinetic 
theory, the characteristic lifetime (e-folding time; x) for a 3 nm particle is: 

X = J]L ( 4 ) 

acA 

The factor y accounts for diffusional resistance to mass transfer that must be taken into account 
for preexisting particles that are on the order of or larger than the air mean free path. For the size 
distr ibutions that we measured in this study, we found that y=1.3±0.2. For the purpose of this 
discussion it is adequate, to assume a fixed value of - 1.3. The average thermal speed, c, for a 
3 nm hydrated H 2 SO 4 particle of density 1.4 g cm‘3 is - 2 x 1 (P cm s * 1 ( assuming T=10°C, 
RH=32%). The preexisting aerosol surface area concentration is A. If we assume the ultrafine 
particles are always captured when they collide with the larger preexisting particles (i.e., a=l), 
the characteristic lifetime of a 3 nm particle at Idaho Hill was — 1-5 hour. This short life span 
likely accounts for the observed decay in numbers of ultrafine particles by afternoon. 

A Representative Day: The H 2 SO 4 * ultrafine particle, and particle surface area 
concentrations on 21 September 1993 are shown in Fig. 3. Like most days at Idaho Hill, H 2 SO 4 



levels began to increase at - 615 MST. The ultrafine particle concentration did not begin to rise 

# 

above nighttime levels until -730 MST, about 1.25 hours later. During the day, the H 2 SO 4 and 
ultra finp. aerosol concentrations tracked well; changes in ultrafine concentrations followed 
changes in H 2 SO 4 concentrations. This general correlation between these species was observed 
on all but one of the sampling days at Idaho Hill. It was also commonly observed at Mauna Loa. 
The measurements suggest that at both sites, H 2 SO 4 was a precursor of the ultrafine particles. 

The measured fine aerosol size distributions for 21 September are shown in Fig. 4. The 
data have been summarized by plotting the average distribution for three periods for which the 
distributions were fairly steady. From Fig. 4a, the largest changes in the number distributions 
were for particles with diameters less than approximately 50 nm. The concentrations of larger 
particles were fairly steady. Note that the concentration of the smallest particles plotted, -20 nm, 
did not significantly increase until approximately 1000 MST. Recall the ultrafine concentration 
(3 nm particles) increased much earlier, at - 715 MST (see Fig. 3). Again, this delay may be the 
time required for growth from -3 nm to -20 nm diameter. 

Another noteworthy feature of the particle size distributions is that at some point during 
the day they had a slight bimodal shape. The minimum between the modes was at -60 nm. This 
bimodal distribution for fine aerosols is a common feature of marine boundary layer aerosols 
where it is understood to result from aerosol processing by non-precipitating marine clouds 
[Hoppel et aL, 1994b]. It was also a regular feature of the fine particle distributions recorded at 
Mauna Loa [Weber and McMurry, 1995]. Though such bimodality has not been commonly 
reported for continental aerosols, at Idaho Hill it was observed at some time during the day for 
nearly one-half the sampling days. 

Figure 4b shows the fine aerosol surface area distributions corresponding to the three 
average number distributions. Plotted in this manner, the area under the curve is the calculated 
aerosol surface area concentration. 

In summary, our data suggests that ultrafine particles were produced by recent nucleation 
and that H 2 SO 4 was a precursor species. The ultrafine particle concentrations varied diumally 



due to the photochemical production of the precursor species and their short life spans. In the 
next section, the influence of parameters expected to affect the nucleation rate are studied by 
constructing scatter plots of ultrafine concentrations versus each parameter. 

Correlations Between Ultrafine Concentrations and Measured Parameters 

For all data free of local contamination and collected under downslope flow, the 
measured ultrafine concentrations are shown in Fig. 5 as a function of the preexisting aerosol 
surface area, relative humidity, and H 2 SO 4 concentration. These plots illustrate the recorded 
range of each parameter and give some indication of their influence on the ultrafine aerosol 
concentration and thus the extent of new particle formation. A point of particular interest is the 
influence of these parameters on the maximum ultrafine concentration since this may provide 
insight about how individual parameters influenced nucleation. In all plots, an envelope roughly 
defining the maximum ultrafine concentration is indicated by a line. Similar scatter plots for 
Mauna Loa data can be found in Weber et aL [1995a]. 

Ultrafine Particle versus Aerosol Surface Area Concentrations: The effect of the 
aerosol surface area on ultrafine concentrations is shown in Fig.5. We found that generally lower 
ultrafine concentrations were recorded during periods of high aerosol surface area 
concentrations. The envelope defining the maximum measured ultrafine concentration decreases 
with increasing surface area. Note that "high" ultrafine aerosol concentrations were never 
dete cte d at Idaho Hill during periods when the aerosol surface area was greater than 
-80 |im 2 cm* 3 . Significant numbers of ultrafine particles were, however, recorded at surface 
areas up to this value. The observed decrease in ultrafine particle concentrations with increasing 
aerosol surface area is consistent with expectations. Concentrations of the nucleating species 
decrease due to heterogeneous condensation as preexisting aerosol concentrations increase; this 
to a reduction in particle production rates. Furthermore, the likelihood that a freshly 
formed nucleus grows to a detectable size before it is scavenged by the preexisting aerosol also 
decrea ses with increasing preexisting aerosol concentrations. 
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Scatter plots of Mauna Loa data {Weber etaL[ 1995a]) show a similar envelope with the 
sa rnft correlation between ultrafine and aerosol surface area concentrations, although the levels of 
both ultrafine particle and aerosol surface area concentrations were much lower at Mauna Loa. 
For example, high ultrafine concentrations were never observed at Mauna Loa during periods 
when the aerosol surface area was larger than -30 jim^ cm'3. Different production rates of 
precursor species (i.e., differing SO2 concentrations) at these two sites probably accounts for 
these observed differences. 

Ultrafine Particle Concentration versus Relative Humidity: From Fig. 5b, the highest 
ultra fine concentrations at Idaho Hill were observed at low relative humidities. Evidence of 
significant particle formation, indicated by high ultrafine particle concentrations, was observecTat 
relative humidities as low as 15 to 30%. At both Idaho Hill and Mauna Loa there is no evidence 
t hat increases in relative humidity led to enhanced particle production, as would be expected if 
the ultra fine particles were formed by binary nucleation of H 2 SO 4 and H 2 O. 

The observed negative correlation between ultrafine particle concentrations and relative 
humidity may be due to the observed positive correlation (correlation coefficient of +0.13) 
between these parameters. This correlation could result from two differing influences. First, 
because the degree of sulfate particle hydration increases with relative humidity, the swelling of 
pa rtic l es with increasing relative humidity leads to higher aerosol surface areas at higher 
humidities; this effect becomes very significant as relative humidies increase beyond 85-90%, 
which is well in excess of values typically encountered in this study. Secondly, relative humidity 
at this site may have been an indicator of air mass origin; more pristine downslope air from 
higher elevations would also likely be dryer. In either case, the negative correlation between 
ultrafine particle concentration and relative humidity likely does not reflect the influence of 
water vapor concentration on atmospheric particle formation rates. 

The observation of particle formation predominantly in dry downslope air (i.e., air 
originating from higher elevations) is similar to the marine observations of Covert et aL [1995]. 
They found evidence for new particle formation only in regions of subsiding air masses. In both 



14 


cases, the lower aerosol surface area concentrations associated with dry air from aloft may have 
been the prerequisite for significant particle production in these regions. 

Ultrafine Particle versus Sulfuric Acid Concentrations: In Fig. 5c, the ultrafine 
particle concentration is plotted with respect to the H2SO4 relative acidity. Sulfunc acid relative 
acidity is the measured H2SO4 partial pressure divided by the saturation vapor pressure of pure 
H2SO4 at the measurement temperature [Ayers et aL , 1980]. Plotted in this way, the major 
influence of temperature on the tendency of H2SO4 and H2O to nucleate is taken into account 
A second abscissa showing H2SO4 concentrations is also shown. In principle, of course, 
concentration and relative acidity cannot be shown on the same graph if temperatures varied 
significantly. In this study the temperature did not vary widely; we used a characteristic value of 
10 *C when evaluating concentrations for this plot From Fig. 5c, it is evident that the measured 
ultrafine concentration and H2SO4 relative acidity were positively correlated; higher ultrafine 
concentrations were, detected at higher H2SO4 concentrations. These results are consistent with 
expectations and our earlier observations [Weber et al. 1995a] that H2SO4 is a primary precursor 
species of the ultrafine particles generated by homogeneous nucleation. 

At both sites, the H2SO4 concentrations typically ranged from -10^ to -lO? cm'3 and 
significant ultrafine concentrations were recorded at unexpectedly low acidities (or H2SO4 
concentrations). Figure 5c shows that at Idaho Hill, high ultrafine concentrations were detected 
at TM' iHitiV-s down to 10‘6 which is much too low for nucleation by H2SO4 and H2O. Another 
major inconsistency with binary nucleation is found when comparing measurements at Mauna 
Loa and Idaho Hill. The maximum ultrafine particle concentrations at a given H2SO4 relative 
acidity were about an order of magnitude higher at Idaho Hill than at Mauna Loa. This suggests 
that species other than H2SO4 and H2O (e.g., ammonia), may also participate in particle 
formation. A conceptual framework for such a nucleation mechanism is discussed by Weber and 
coworkers [1995b; 1995c]. 
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Discussion 

Ultrafine Particle Growth Rates 

In this section, the data are used to estimate average growth rates of sub-3 nm particles. 
These rates are compared to growth by condensation of only hydrated H 2 SO 4 vapor calculated 
using Eq. (1) and the measured H 2 SO 4 concentration and relative humidity. Two approaches are 
used to estimate the growth rate of ultrafine particles. First, growth rates are inferred from the 
observed time lag between the rise in H 2 SO 4 concentration and the rise in ultrafine particle 
concentrations after sunrise. In the second approach, the observed effects of cluster scavenging 
by the preexisting aerosol surface area on ultrafine particle concentrations were used to estimate 
the ultrafine particle growth rate. 

Growth Rates Inferred from the Delay Between H 2 SO 4 and Ultrafine Particle 
Concentrations At Sunrise: As shown in Fig. 3, the H 2 SO 4 concentration began to increase at 
-615 MST on 21 September 1993 whereas the ultrafine concentration began to rise at -730 
MST. If we assume this delay is due solely to the time required for a H 2 SO 4 vapor molecule of 
rfiamrtw -1 nm to reach our lower detection limit of -3 nm, then the observed average growth 
rat e was 2 nm in 1.25 hours, or -1.6 nm h _1 . From the measurements of H 2 SO 4 and H 2 O 
concentrations between 615 and 730 MST, Eq. (1) predicts an average growth rate from 
condensing H 2 SO 4 vapor and its associated water of -0.2 nm h* 1 which is a factor of 8 below 
the observed growth rate. On this morning, if growth was solely due to condensation of H 2 SO 4 
vapor, the ultrafine particle concentration would not have begun to increase until about 1045 
MST. 

Recall that later in the morning of 21 September the concentrations of 20 nm particles did 
not increase until -1000 MST (Fig. 4a). The number of ultrafine particles (-3 nm) had increased 
about 2.5 hours earlier at 730 MST. For 3 nm particles to reach 20 nm in 2.5 hours requires an 
average growth rate of -7 nm h' 1 . During this time the average growth rate by Eq. (1) based on 
measured H 2 SO 4 concentrations was -0.6 nm h* 1 ; the observed growth rate was roughly 1 1 
times higher. 



Observed growth rates for other days are compared in Table 1 with values calculated 
from Eq. (1). Results are shown for all days for which we had both H 2 SO 4 and ultrafine aerosol 
concentrations in the morning. Note that growth rates of sub-3 nm particles were consistendy 8 
to 13 times higher than growth by condensation of only hydrated H 2 SO 4 vapor. 

Growth Rates Inferred From The Effects of Cluster Scavenging: Insights about 
particle growth rates can also be obtained by investigating the relationship between ultrafine 
particle concentrations and aerosol surface area. The probability that a freshly nucleated particle 
will grow to our minimum detectable size (-3 nm) decreases with increasing aerosol surface 
area. This probability also decreases with decreasing particle growth rates: lower growth rates 
lead to longer growth times during which losses to preexisting particles can occur. 

The time-dependent concentration, N, of a population of uniformly-sized particles that is 
growing by gas-to-particle conversion and is being scavenged by preexisting particles can 
approximately be expressed as: 

dN _ Nc A (5) 

dt ” 4 y 


where c is the mean thermal speed of the particles and A is the surface area concentration of the 
preexisting particles. As was discussed above, we assume 7 = 1 - 3 . From kinetic theory, the mean 
thermal speed for spherical particles of diameter D p and density p is: 


c(D P ) = 


48kT 

^pDp j 


( 6 ) 


where k is Boltzmann’s constant, and T is absolute temperature. As was discussed above 
(Eq. (1)), condensational growth rates are independent of size if evaporation is negligible relative 



to con dens ation Assuming that growth rates, dDp/dt are also independent of time, the time 
dependent size is: 


D p (t) = D pl+ ^-. W 

where Dpi is size of the freshly nucleated particle. Substituting (6) and (7) in (5) and solving for 
N(Dp 2 (t» leads to the following result for the probability, P, that a particle will grow from Dpi 
to Dp2 before it is scavenged by the preexisting aerosol: 



( 8 ) 


Freshly nucleated particles probably consist of a few molecules. We will assume that D p i = 

1 nm; our discussion is not particularly sensitive to the assumed value. Furthermore, we assume 
that Dp 2 = 3 nm, our minimum detectable size. It follows that: 


p N (°p2) 
N(D pl ) 



-6.8xl0“ 2 


* 

A 

dDp/dt * 


(9) 


where A is in pm^m* 3 , dDp/dt is in nm h* 1 , and a density of 1.4 g cm* 3 was assumed. If growth 
is due solely to condensation of H 2 SO 4 vapor, then from (1) and (9): 


_ N(D 2 ) 

P c— = exp 

N(D pl ) 


-9.2 xlO 5 - 


n h 2 so 4 J 


( 10 ) 


where Nh 2 so 4 “ die concentration of gas phase sulfuric acid in molecules cm* 3 . An effective 
molecular weight of the condensing sulfuric acid and its associated water, etc., of 200 was 


assumed. 



Figure 6 shows the relationship between the peak ultrafine particle concentration 
m eas ured on each day at Idaho Hill and the preexisting aerosol surface area concentration. A. 
B ased on the hypothesis that particles are produced by a multicomponent nucleation process that 
involves sulfuric acid, data are categorized by the peak sulfuric acid concentration that was 
measured at or prior to (typically within 1 to 2 hours) the time that the peak in ultrafine 
concentrations was observed. We confine our attention to the peak daily ultrafine concentrations 
dnrv» numerical calc ulati ons justify assuming steady state cluster size distributions that are 
dependent on the concentration of the nucleating vapor at this point (e.g., (McMurry, 1983; Rao 
and McMurry, 1989)). Note that there is a general trend towards decreasing peak ultrafine 
concentrations with increasing aerosol surface areas and with decreasing sulfuric acid 
concentrations, as expected. 

Also shown in Figure 6 are two lines obtained from Equation 10 assuming 
n H 2 S0 4 = 3*10 6 and lOxlO 6 cm -3 . Because we do not know how many particles are produced 

at a given sulf uric acid concentration, we have arbitrarily set N(0) = 10 3 cm -3 for both of these 
curves; our interest here is in comparing the slopes of these curves with trends in the data. 
Although there is too much scatter in the data to obtain accurate values for the logN versus A 
slopes, it appears that they are ~5 to 10 times smaller than predicted by Equation 10. This 
suggests that growth rates are 5 to 10 times higher than can be explained by H 2 SO 4 
condensation. This is consistent with the results in Table 1 that were obtained independendy, 
and further supports our finding that a species in addidon to H 2 SO 4 is contributing to growth of 
freshly nucleated ultrafine panicles. Our work ( Marti et al, 1995) showed no relationship 
between concentrations of ultrafine particles and of alpha or beta pinene (or estimates of their 
reaction products). This suggests that, at this site, terpenes played a minor role in new particle 
formation compared to sulfates. Terpenes, however, were found to be correlated with the aerosol 
surface area and volume concentrations. This may indicate that these species contributed to 
particle growth. If the high growth rates were from heterogeneous condensation of organics, a 
significant fraction of the ultrafine particles would be organic. 



Predicted H2SO4 Concentrations 

If we assume that the only source of sulfuric acid is the OH, SO2 reaction, and that the 
primary removal mechanism is condensation on preexisting particles with an accommodation 
coefficient of 1.0, the steady state sulfuric acid concentration is: 

[H 2 S0 4 ]= ™°* 3 (11) 

F 

where k is the rate constant for the SO2, OH reaction and F • [H2SO4] equals the rate of sulfuric 
acid transport to preexisting particles. We take k = 1.1x10*12 cm^s'l [DeMore et aL, 1992; 
Gleason et aL, 1987; Wang et aL, 1988], The loss factor, F, was obtained by integrating the 
transition regime expression of Fuchs and Sutugin [1970] across the measured aerosol size 
distribution. The steady state assumption is justified since the characteristic time required to 
establish steady state is relatively short (-3 min). Sulfuric acid concentrations predicted with 
Fq (11) are compared with measured values for 24 July 1992 (Mauna Loa) and 21 September 
1993 (Idaho Hill) in Fig. 7. These two days were selected because, unlike most days, the 
recorded SO2 concentrations were sufficiently above the instrument’s lower detection limit to 
provide confidence in the measured values. Median daytime H2SO4 production rates were 6x10^ 
and 3x10* craV 1 at Mauna Loa (7/24) and Idaho Hill (9/21) respectively. Median fractional 
H2SO4 loss rates on these days were were 5xl(H and 6xl0* 3 s* 1 . Measured and calculated 
steady state concentrations are in very good agreement 

The loss Factor, F, in Eq. (1 1) is an implicit function of the H2SO4 mass accommodation 
coefficient on aerosol particle surfaces. In our analysis we have assumed that the 
accommodation coefficient equals 1.0. Previous researchers have used values for mass 
accommodation coefficients ranging from about 0.3 to 0.04 [Hegg et aL, 1990; Roes et al, 1992 
and 1995]. These values would lead to increases in predicted steady-state sulfuric acid 



concentrations by factors of -3 to - 20 . We conclude that given the good agreement between 
measured and calculated values, the assumed value of 1.0 is in the right range. 

In deriving Equation ( 1 1 ) we also implicitly assumed that the equilibrium vapor pressure 
of sulfuric acid above the aerosol particles is small relative to the measured sulfuric acid 
concentrations. If this had not been the case, then the measured sulfuric acid concentrations 
would have systematically exceeded the calculated values. Again, there is no evidence to 
suggest that reevaporation was a significant source of sulfuric acid vapor on these days. 

Conclusions 

Measurements at a remote continental site indicate that elevated concentrations of 
ultrafine particles (nominally 3-4 nm diameter) resulted from recent new particle formation. The 
data points to H2SO4 as a precursor vapor of these newly formed particles. 

The measured ultrafine particle concentration had a consistent diurnal pattern. There is 
evidence that this was due to the photochemical production of the precursor species (i.e., H2SO4 
and possibly others) and the relatively short lifetime (- 1.5 hour) of the ultrafine particles. 
Although the H2SO4 concentration was observed to increase just after sunrise, elevated ultrafine 
particle concentrations were delayed by approximately 1 hour. Growth rates of ultrafine particles 
were estimated from this delay and found to be roughly 8 to 12 times higher than growth by 
condensation of hydrated H2SO4 vapor. By a completely different approach, ultrafine particle 
growth rates were estimated from the influence of cluster scavenging on ultrfine particle 
concentrations. With this method, we estimate that growth rates of ultrafine particles were 
roughly 5 to 10 times higher than growth by H2SO4 vapor. 

New particle formation, indicated by high ultrafine concentrations, was common at Idaho 
Hill in downslope air and was observed during periods of low relative humidity (median RH was 
32 %). Significant ultrafine concentrations were measured at H2SO4 relative acidities as low as 
10*6. These are water vapor and acidities for which classical H2SO4-H2O nucleadon theory 
would predict practically no nucleation. 



At both Mauna Loa and Idaho Hill, the two major factors influencing particle formation 
were found to be the H 2 SO 4 and preexisting aerosol surface area concentration. Sulfuric acid 
appeared to be a primary precursor species of the ultrafine particles and the highest ultraflne 
concentrations were recorded when the aerosol surface area was low (relative to typical surface 
area concentrations for that particular site). For a selected day at each site, the steady state 
H 2 SO 4 concentration, calculated from the balance between H 2 SO 4 photochemical production 
and H 2 SO 4 scavenging by aerosol, agreed remarkably well with the measured values. 

Compared to Mauna Loa, new particle formation was more vigorous at Idaho Hill despite 
similar H 2 SO 4 relative acidities at both locations and much lower relative humidities at Idaho 
Hill We speculate that additional species, such as ammonia, may have also participated in new 
particle production and that higher levels of these species at the continental site enhance 
nucleation there. Participation of ammonia, however, is not expected to significantly enhance 
particle growth rates, although other species, such as organics, may. These species may not have 
been involved in new particle formation, but could significantly enhance the growth rates of the 
newly formed particles and account for the high ultrafine particle growth rates observed at Idaho 
Hill. 
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Table 1 : Comparison of estimated observed growth rates of ultrafine particles and calculated 
growth rates based on measured H2SO4 and H2O vapor concentrations assuming that 
growth is by condensation of hydrated H2SO4 vapor. The ratio of the observed 
growth rates to the calculated rates is shown in the right most column. 


Day of Sept. 
1993 

Observed 

Growth, 

(nra/h) 

Calculated 
Growth, . 
(nra/h) 

Ratio, 

ObsJCalc. 

6 

1 

0.1 

10 

10 

1.3 

0.2 

9 

11 

1.1 

0.1 

12 

12 

1.3 

0.1 

12 

15 

2 

0.1 

10 

21 

1.6 

0.2 

8 

23 

1.3 

0.1 

12 

24 

1.1 

0.09 

12 

26 

0.5 

0.04 

13 



Figure Captions 


Figure 1 Total aerosol number concentrations (diameters larger than 3 nm) (a), and ultrafine 

aerosol concentrations (diameters between nominally 3 nm and 4 nm) (b), as a function 
of wind direction. The plot shows that high ultrafine aerosol concentrations were 
observed predominantly in air from the western sector (downslope). 

Figure 2 Recorded daily concentrations of OH (a), SO 2 (b), H 2 SO 4 (c), and ultrafine particles 
(d). For all plots, only data during periods of downslope air flow is shown. Also 
plotted is the UV solar intensity on a clear day (26 September 1993). Both OH and 
H 2 SO 4 concentrations often increased from nighttime levels at sunrise. Ultrafine 
aerosol concentrations, however, never began to rise until over an hour later. This 
delay is used to estimate the growth rate of ultrafine particles. 

Figure 3 Measured H 2 SO 4 (a), ultrafine aerosol (b), and aerosol surface area (c) concentrations 
on 21 September 1993. On this day the air flow was downslope from 000 through to 
1700 MST. 

Figure 4 Time average fine aerosol number and surface area distributions for periods 600 - 
1000, 1000 - 1500, and 1500 - 1800 MST on 21 September 1993. Increased 
concentrations of the smallest particles did not begin until - 1000 MST; concentrations 
of the higher mode were fairly steady. The fine aerosol also appears to be bimodal; 
similar (albeit more strikingly bimodal) distributions have been reported for remote 
marine aerosols. 

Figure 5 S catter plots of the measured ultrafine concentration as a function of the aerosol 

surface area (a), the relative humidity (b), and the H 2 SO 4 relative acidity (c). In plot 
(c), the corresponding H 2 SO 4 concentration is shown assuming a temperature of 10°C. 
All plots contain data for periods of only downslope air flow. 
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Figure 6 Maximum ultrafine particle concentration observed on individual days versus 

preexisting aerosol surface area concentrations on those days. The data are categorized 
by the sulfuric acid concentrations that were measured about one hour prior to the peak 
ultrafine concentrations on the hypothesis that this was responsible for new particle 
production. The slopes of the solid lines (from Eq. (10)) show the expected sensitivity 
of ultrafine concentrations to aerosol surface area if growth rates are limited by H2SO4 
compensation. 

Figure 7 Comparison of the measured H2SO4 concentration and the predicted steady state 
concentration for 24 July 1992 at the Mauna Loa Observatory, Hawaii, and 21 
September 1993 at Idaho Hill, Colorado. The predicted steady state H2SO4 
concentration is determined from the balance between its production by the SO2 - OH 


reaction and loss by aerosol scavenging. 
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ABSTRACT 


Measured production rates of tropospheric ultrafine particles (-3 run diameter) are 
reported for the first time and are shown to be orders of magnitude greater than 
nucleation rates predicted by the binary theory of homogeneous nucleation for sulfuric 
acid and water. Furthermore, the functional dependence of observed particle formation 
rates on sulfuric acid vapor concentrations is much weaker than predicted by binary 
theory. We present arguments to show that these discrepancies might be due to the 
participation of a species such as ammonia which could stabilize subcritical clusters, 
thereby enhancing nucleation rates. The data suggest that atmospheric nucleation may 
occur by a collision-limited process, rather than by a condensation/evaporation controlled 
process as is assumed in the classical theory. 


INTRODUCTION 


New particle formation by nucleation of gas phase species significantly influences 
the size and numbers of tropospheric aerosols. These aerosols can affect the earth's 
radiation budget directly by scattering solar radiation (Charlson et aL, 1992; Kiehl and 
Briegleb, 1993) or indirectly by serving as cloud condensation nuclei (Twomey, 1977). 
Since in both cases the magnitude of the effect is sensitive to particle size and number 
density, accurate representation of nucleation is essential for accurate modeling of 
aerosol-climate coupling. 

Atmospheric particle nucleation has been studied in the laboratory and simulated 
by numerous modeling studies, but in situ studies are more limited. Several recent field 
measurements have demonstrated without question that the formation of new particles by 
homogeneous nucleation occurs in the atmosphere. Table 1 lists work published since 
1990 that presents evidence of new particle formation in a variety of environments. 

Much attention has been paid to nucleation in the marine boundary layer, and to a lesser 
extent, the marine free troposphere. Particle formation near marine clouds has received 
particular scrutiny. Nucleation over remote continental areas remains relatively 
unexplored. 

Recent advances in instrumentation for detecting ultrafine particles has provided 
new insights into atmospheric nucleation phenomena. In particular, the ultrafine 
condensation particle counter (UCPC), developed in our laboratory (Stolzenburg and 
McMurry, 1991) and commercialized by TSI (St. Paul, MN) as the UCPC Model 3025, 
has played an important role in this work. With a 50% detection limit of less than 3 nm, 
the UCPC allows new particles to be counted sooner after nucleation than was previously 
possible. 

Presently, despite little evidence for its applicability, most modeling studies of 
atmospheric aerosol nucleation and growth (e.g., Kreidenweis et aL, 1991; Raes et aL, 
1992; Raes and van Dingenen, 1992; Hegg, 1993; Lin et al., 1993; Easter and Peters, 

1994; Russell et al., 1994; Raes, 1995) use the classical binary theory of homogeneous 
nucleation for sulfuric acid (H2SO4) and water (H2O) (Jaecker-Voirol and Mirabel, 

1989). According to classical theory, new particles are bom when molecular clusters 
grow past a stable (critical) size (Shugard et al. 1974). In binary H2SO4/H2O nucleation, 
cluster growth rates are determined by the rates at which H2SO4 and H2O are transported 
to and evaporate from subcritical clusters that contain only H2SO4 and H2O. Because the 
atmosphere contains a complex spectrum of species that could participate in the 
nucleation process, there is no a priori reason to expect the H2SO4/H2O theory to be 



generally valid. 

In this paper we show that actual atmospheric rates of new particle formation 
measured at Mauna Loa, HI and Idaho Hill, CO are far faster than predicted by the theory 
for binary nucleation of H 2 O and H 2 SO 4 . Furthermore, the observed dependence of 
particle production rates on H 2 SO 4 concentrations is inconsistent with classical 
theoretical expectations, but is more nearly consistent with a collision-controlled process. 
We provide evidence to show that NH 3 is likely to participate in nucleation, and may be 
responsible for the observed high rates. 

When Peter McMurry began his doctoral research under the supervision of 
Sheldon Friedlander at Caltech in 1972, Friedlander suggested that he work on the 
problem of nucleation in chemically reacting systems. This problem has been an 
enduring theme of McMurry's research since that time. Chemical nucleation is important 
in the atmosphere, in semiconductor processing equipment, and in aerosol reactors used 
for material synthesis. While a great deal of elegant work has been done to evaluate the 
classical liquid droplet model for nucleation, it is well known that there are significant 
discrepancies between theory and experiment remain. Furthermore, the multicomponent 
thermodynamic data needed to apply this theory to aerosol systems of practical 
importance is often unavailable. This is because vapor pressures for species that nucleate 
under many practical circumstances are exceedingly low and difficult to measure. In 
addition, nucleation frequendy occurs under circumstances for which the classical theory 
is invalid. For example, if the vapor pressures of chemically-produced condensable 
species are sufficiendy low, then the nucleation process will be collision-controlled rather 
than condensation/evaporation controlled as is assumed in the classical theory. 
Alternatively, the classical theory is invalid when particle growth is dominated by surface 
reactions (as likely occurs in semiconductor processing equipment) rather than by 
condensation. Validated approaches for modeling such phenomena are generally not 
available. The present paper summarizes our most recent thinking on atmospheric 
nucleation. Our ultimate objective in this research is to obtain data that can be used to 
develop and verify theories for nucleation in the atmosphere. 

EXPERIMENT 

In this paper we estimate particle formation rates from measurements made at a 
remote marine site (Mauna Loa Observatory, Hawaii 1992) and a remote continental site 
(Idaho Hill, Colorado, 1993). These estimates were made possible by the recent 
development of instrumentation for measuring concentrations of gas phase H 2 SO 4 and 



3 to 4 nra ultrafine particles. The H2SO4 was measured using an atmospheric pressure- 
selected ion chemical ionization mass spectrometric technique (Eisele and Tanner, 1993). 
The ultrafine particles were measured using the ultrafine condensation particle counter 
(UCPC) described by (Stolzenburg and McMurry, 1991) in the pulse height analysis 
(PHA) mode. 

The PHA technique has been described previously (Brockmann, 1981; 
Stolzenburg, 1988; Ahn and Liu, 1990; Wiedensohler etaL, 1994); the reader is referred 
to these previous publications for details. In summary, particles larger than -15 nm are 
activated (i.e., begin to grow by condensation) near the entrance to the UCPC condenser 
and all grow to about the same final droplet size. However, curvature effects become 
significant for particles smaller than this. Therefore, such "ultrafine" particles must travel 
further into the UCPC condenser before saturation ratios are sufficiently high to initiate 
their condensational growth. Because growth times decrease with initial size, the final 
droplet size at the exit from the condenser also decreases with initial size. This final 
droplet size can be measured with the photo detector that is used to count individual 
particles exiting from the condenser. Thus, by measuring the UCPC pulse height 
distributions it is possible to obtain quantitative information on ultrafine particle size 
distributions. 

For the data reported in this paper, counting intervals were 15 minutes (Mauna 
Loa) and 2.5 minutes (Idaho Hill). We confine our attention to a range of pulse heights 
that excludes all particles exceeding 4.5 nm. Because the minimum detectable size for 
this instrument is -2.7 nm, we assume that the mean size of the counted particles is 
-3.5 nm. We also assume that all particles in the 3 to 4 nm diameter range are counted, 
although we know from measured instrument PHA response functions that a fraction of 
the particles in this range will produce pulse heights that exceed our lower limit Because 
of these factors, our reported particle formation rates are somewhat smaller than actual 
values. Uncertainties introduced by these approximations , however, are likely to be 
small compared to the variability in the data. A more complete discussion of the PHA 
technique that was used for these measurements is discussed by (Weber, 1995). 

The PHA method offers considerable advantages over methods mentioned in 
Table 1 for measuring size distributions and concentrations of ultrafine atmospheric 
particles. Techniques that rely on calculating the difference between concentrations 
measured with two CPCs having differing lower detection limits (e.g. Nuitrafine=N3025- 
N3760) are severely limited by counting-statistical uncertainties and flow uncertainties. 
For example, if the concentration of particles larger than 15 nm is -200 cm* 3 (a typical 
value in the remote troposphere) and if a time resolution of 5 minutes is required, then 



counting statistical uncertainties limit the minimum ultrafine (3-15 nra) concentration that 
can be measured with 15% accuracy to ~>10 cm'3. This value will increase when the 
flow uncertainties of the two instruments are taken into account In contrast using the 
PHA technique under the same sampling condidons, ultrafine concentrations as small as 
0.1 cm~3 can be measured to within -10-15%. Furthermore our measurements are only 
weakly dependent on flow uncertainties, and provide information on size as well as 
concentration. Techniques that employ differential mobility analysis (Le., DMPS or 
SMPS systems) can provide excellent size resolution, but are limited by counting rates: 
only a very small fraction of the smallest particles are charged, and only charged particles 
are classified and counted with this technique. Thus, far more time is require to measure 
size distributions. Diffusion batteries require time to scan through a series of sampling 
ports, and also require a steady size distribution for accurate measurements. In summary, 
the PHA technique provides information on 3 to 12 nm size distributions with better time 
resolution and sensitivity than is possible with other methods. 


THEORY 


In this section the approach that is used to estimate rates of 3-4 nm particle 
formation from our field measurements is discussed. The rate at which particles grow 
past the minimum detectable size (-3 nm diameter) by vapor condensation is: 


J(3 nm) = 


dN 

dD p 


13.5 nm 


dDp 

dt 


3.5 nm 


[ 1 ] 


where dN/dDp is the measured aerosol size distribution function at 3.5 nm and dDp/dt is 
the diameter growth rate for 3.5 nm particles. The average size distribution function in 
the vicinity of 3.5 nm, n, can be estimated from our ultrafine measurements by: 


dN 

dD p 


13.5 nm 


= n(3.5 nm) 


measured concentration in 3 to 4 nm diameter range 
4 nm - 3 nm 


•[ 2 ] 


We have previously shown that the appearance of ultrafine particles is correlated 
with H2SO4 vapor concentrations, suggesting that H2SO4 participates in new particle 
production (Weber et al., 1995). If condensation of H2SO4 vapor is much faster than 
evaporation (a reasonable assumption in our case where particle size greatly exceeds the 



critical size for binary H2SO4/H2O homogeneous nucleation), and if the rate limiting step 
for particle growth is uptake of H2SO4, then from kinetic theory (e.g. T Present, 1958 ): 


dv 

dt 


collisions total volume added 


area time 


collision 


particle surface area 


[H 2 SQ 4 ]c 


jcD? dD, 


2 

= «h 2 so, • • 4 - V 1*D„ — 


dt 


Therefore: 


[ 3 ] 


[ 4 ] 

where [H2SO4] is the molecular gas phase concentration of H2SO4, a H ,so 4 is the H2SO4 

mass accommodation coefficient, c is the mean thermal speed of the condensing H2SO4 
species, and v l is the amount that the particle volume increases upon the addition of a 
single H2SO4 molecule. This volume increment includes species that are transported to 
the particle with the H2SO4 as well as species that are absorbed afterwards to maintain 
phase equilibrium, such as H2O and possibly NH3. Actual diameter growth rates will 
exceed the value given by Eq. [ 4 ] if additional species that are not associated with H2SO4 

are involved. Assuming that every collision between 3 nm particles and the hydrated 
H2SO4 molecule is effective (a HjS0< = 1 ), the "measured” rate of new particle formation 

is: 


Jmeasured^n-l^i^^c. [ 5 ] 

Note that rates calculated by Eq. [ 5 ] are based on instantaneous point data, and are 
therefore not affected by advection. 

RESULTS AND DISCUSSION 

Figure 1 shows particle formation rates calculated using Eq. [ 5 ] versus measured 
H2S0 4 concentration for all data obtained at Mauna Loa and Idaho Hill that were judged 
free of local anthropogenic influences. In both cases the data are bounded by two lines. 
The maximum measured particle formation rates vary approximately with the square of 
the H2SO4 concentration (line of slope -2 in Figure 1 ). No significance regarding 
nucleation can be drawn from the apparent lower boundary (line of slope 1) which is an 



artifact of the measurement procedure. This lower boundary corresponds to 
measurements in which only one particle was detected by the UCPC during the 
measurement interval. Sulfuric acid concentrations for cases when no ultrafine particles 
were detected are given at J = 10*10 cm* 2 s* 2 . 

Predicted rates of homogeneous nucleation according to the classical hydrate 
theory for binary (H2SO4/H2O) nucleation are also shown in Figure 1 . The measured 

particle formation rates significantly exceed the theoretical rates. While classical theory 
predicts practically no nucleation for [H2SO4] < 10 ^ cra'l (at 50 % RH), we observed new 
particle formation at H2SO4 concentrations down to 10 4 cra*l. Because the critical size 
for homogeneous nucleation (~ 1.5 nm diameter) (Mirabel and Katz, 1974 ) is significantly 
below our minimum detectable size (-3 nm), the predicted rates of homogeneous 
nucleation are not directly comparable with our reported rates of new particle formation. 
In particular, a fraction of the newly-formed particles is scavenged by preexisting 
particles before they grow to 3 nm. Therefore, if binary theory were used to model 
aerosol production, predicted rates of 3 nm particle production would be lower than the 
theoretical rates of binary nucleation shown on Figure 1 . Thus scavenging by preexisting 
particles will further increase discrepancies between our measured rates, which include 
scavenging, and the theoretically-predicted values which do not 

Although the experimental values in Figure 1 were obtained assuming a H2SO4 
mass accommodation coefficient of 1.0, the discrepancy between theory and experiment 
would not be resolved if the actual value were smaller than this. The accommodation 
coefficient would not affect the functional dependence of our reported rates on H2SO4 
vapor concentrations. Furthermore, nucleation rates from binary theory would also need 
to be decreased to account for nonaccommodation, thus not resolving the discrepancy 
between theory and measurements. 

The dependence of new particle formation rates measured in Colorado on aerosol 
surface area is shown in Figure 2 . Note that higher aerosol surface areas led to reduced 
rates of particle production. This is consistent with expectations, since clusters formed by 
nucleation can be scavenged by preexisting particles before they grow to our minimum 
detectable size. Thus, aerosol surface area accounts for some but not all of the scatter in 
the Colorado particle formation rate results. A similar dependence on surface area was 
not observed in Hawaii. This may be because of the higher measurement uncertainties 
associated with the lower concentrations in that environment 

The observation that the maximum measured new particle formation rates vary in 
proportion to [H2SO4] 2 is intriguing. This suggests that the nucleation process may be 
collision-controlled rather than evaporation/condensation controlled, as is assumed in 



classical nucleation theory. McMurry (1980) has discussed the problem of aerosol 
formation in a reaction system when nucleation is collision-controlled, and Rao and 
McMurry (1989) have discussed the transition between the collision-controlled and 
condensation/evaporation-controlled limiting regimes Assuming a steady-state cluster 
distribution, the upper limit for the new particle formation rate is the rate at which H2SO4 
collides with itself, J3[H2SC>4] 2 , where the collision frequency function (3 is about 3x1 (H^ 
cm 3 s* 1 (McMurry and Friedlander, 1978; McMurry, 1983a). Observed maximum new 
particle formation rates at Hawaii and Colorado were, respectively, approximately six and 
five orders of magnitude below this, (Figure 1). Phenomena that will reduce rates below 
the H2SO4 collision rate include cluster evaporation, nonaccommodation for H2SO4 
collisions, cluster scavenging by preexisting aerosol as clusters grow up to our minimum 
detectable size (-3 nm), or the participation of species other than H2SO4 and H2O in the 
nucleation process. The latter is the only explanation consistent with our data. The 
binary theory, while accounting for the effects of cluster evaporation, produces results 
that are inconsistent with our measurements; including nonaccommodation in the binary 
theory would increase discrepancies with our data. The Colorado data in Figure 2 show 
that cluster scavenging by preexisting particles led to reduced rates of particle production 
at given H2SO4 concentrations during our measurements, but not when new particle 
formation rates were maximum (i.e., along the J~[H2S04] 2 line shown in Figure 1). 

Scenarios that will lead to the observed results can be found if a stabilizing 
species such as NH3 were to react with clusters containing one or more H2SO4 molecules 
and H2O. Such species would decrease rates of H2SO4 evaporation from clusters, 
thereby increasing nucleation rates. For example, nucleation rates could be determined 
by collision rates of a thermodynamically-determined NH3-stabilized fraction, y, of the 
H2SO4 monomer. The resulting steady state rate of particle formation in the absence of 
cluster scavenging by preexisting particles would then be: 

^L-Pr 2 [H 2 so 4 ] 2 [<n 

where (3 is the collision frequency function for the molecular clusters that contain the 
stabilized H2SO4 molecule. Assuming that p - 3xlO* 10 cm 3 /s, it follows from the data in 
Figure 1 that y-0.001 for Hawaii and 0.003 for Colorado. Alternatively, stabilization 
might occur at a later step in cluster growth. In this case only a fraction of the clusters 
produced by H2SO4 collisions would be stabilized and begin growth towards the 
detectable size. Unfortunately, NH3 was not measured during the field studies discussed 



here. Previous measurements have shown, however, that NH3 concentrations in regions 
similar to those at Idaho Hill are roughly an order of magnitude higher than values found 
over the Pacific (Ayers and Gras, 1983; Tanner and Eisele, 1991; Langford and 
Fehsenfeld, 1992). It is interesting to note that the maximum particle formation rates 
shown in Figure 1 are also about an order of magnitude higher in Colorado than in 
Hawaii. 

There is good reason to believe that NH3 may participate in atmospheric 
nucleation (Hoppel, 1975; Scott and Cattell, 1979). Characteristic values for 
concentrations of H2O, NH3, and H2SO4 in the remote troposphere are ~10 17 cm* 3 , 10^ 
to 10 10 cm -3 , and 10 4 to 10 7 cm' 3 , respectively. Because the NH3 and H2O are typically 
greatly in excess of H2SO4, these species collide with molecular clusters at much higher 
rates than does H2SO4. Laboratory experiments have shown that the reaction probability 
for NH3 with aqueous H2SO4 droplets is close to 1.0 (Huntzicker et al., 1980; McMurry 
et al., 1983b); it is likely that NH3 also reacts efficiently with atmospheric H2SO4/H2O 
clusters. Furthermore, we recently completed laboratory experiments which showed that 
vapor pressures of H2SO4 over H20/H2S04/(NH4)2S04 solutions drop by orders of 
magnitude as ammoniura-to-sulfate ratios rise above one (Marti et al., 1995). These 
observations suggest that tropospheric NH3 is probably present in sufficient quantities to 
stabilize clusters, thereby leading to higher nucleation rates than would occur in its 
absence. It is also possible that additional species not considered here play a role in 
nucleation. Our previous work has shown, however, that one candidate species, methane 
sulfonic acid, is not correlated with ultrafine particles at Mauna Loa (Weber et al., 1995). 

Our field measurements have the limitation that they were made at a fixed 
location. This can complicate our interpretations because there is about a 1 hour growth 
period between the time particles are produced by homogeneous nucleation and the time 
they are detected with our instrumentation. Therefore, the H2SO4 concentration present 
during the nucleation event could have been different from the H2SO4 concentrations 
measured at our site. Although this would not affect the local particle formation rates 
calculated with Equation [5], it could affect the inferred functional dependence of particle 
formation rates on H2SO4 concentrations. However, supplemental measurements 
strongly suggest that homogeneous nucleation was occurring at the measurement site. A 
mass spectrometer was used to detect sub-3 nm (2xl0 3 to 10 4 arau) molecular clusters for 
a portion of the measurement period at Mauna Loa (Weber et al., 1995). We found that 
clusters were typically present when ultrafine particles were detected. Furthermore, at 
both sites, the observed delay times between increases in H2SO4 concentrations and the 
appearance of 3 nm particles were often comparable to expected values, based on 



diameter growth rates given by Equation [3]. Therefore, while transport may have 
contributed to scatter in the observed relationships between particle formation rates and 
H 2 SO 4 concentrations shown in Figure 1, the evidence that some of the measurements 
were made during nucleation events and the consistency of the results obtained in two 
distinctly different locations supports our conclusion that the results of Figure 1 were 
determined primarily by local transformations and not by transport 

CONCLUSIONS 

In summary, our data shows that although H 2 SO 4 is involved with the formation 
of new particles in the troposphere, observed rates of new particle formation typically 
exceed rates of homogeneous nucleation predicted by the classical hydrate theory for the 
H 2 SO 4 /H 2 O system. Furthermore, particle formation is observed at lower H 2 SO 4 
concentrations than predicted by binary nucleation theory. We postulate that these 
discrepancies between observation and theory can be attributed to the stabilizing 
influence of species such as NH 3 . Ammonia is known to react at a near collision- 
controlled rate with H 2 SO 4 /H 2 O droplets, and markedly reduces equilibrium H 2 SO 4 
vapor pressures when incorporated into H 2 O/H 2 SO 4 solutions. It is also possible that 
other species not considered here also participate in the nucleation process. Thus, 
because new particle formation plays a key role in regulating aerosol size distributions 
and concentrations of cloud condensation nuclei, species such as NH 3 may significantly 
influence our view of the dynamic interaction between aerosols and climate. Future work 
incorporating these results into atmospheric aerosol models may lead to new insights 
regarding the impact of aerosols on climate. 
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FIGURE CAPTION: 

Figure 1 Measured rates of new particle formation at a remote marine site, Mauna Loa 
Observatory, Hawaii (a), and a remote continental site, Idaho Hill, Colorado, 
located on the eastern face of the Rocky Mountains (b). The median relative 
humidities and temperatures at Mauna Loa and Idaho Hill for the data shown 
were: 48%, 8 °C and 32%, 6 °C respectively. The measured particle formation 
rates are compared to classical hydrate binary (H 2 SO 4 /H 2 O) nucleation theory 
for relative humidities of 100, 50 and 10%, and with the collision rate of 
hydrated (50% RH) H 2 SO 4 molecules, . 

Figure 2. Dependence of averaged particle formation rates (from Figure lb) on 
preexisting aerosol surface area concentrations. 



Table 1 : Recent field measurements of new particle formation . 


Location 

Reference 

Comment 

Arctic marine boundary 
layer 

Covert et al.{ Covert et 
aL, 1995) 


Marine boundary layer 

Weber et al. (Weber, 
McMurry et al., 1995) 

Simultaneous CN and 
H 2 SO 4 measurements 

Marine boundary layer 

Hoppel et aL (Hoppel et 
al., 1994) 

In clear air and near 
stratiform clouds 

Marine free troposphere 

Perry and Hobbs (Perry 
and Hobbs, 1994) 

In clear air, adjacent to 
marine cumulous clouds 

Arctic marine boundary 
layer 

Wiedensohler et aL 
(Wiedensohler et al., 
1994) 

Intercomparison of 4 
ultrafine aerosol 
measurement techniques 

Remote continental 
boundary layer 

Koutsenogii and 
Jaenicke (Koutsenogii 
and Jaenicke, 1994) 


Marine upper 
troposphere 

Clarke (Clarke, 1993) 


Marine boundary layer 

Putaud et al. (Putaud et 
al., 1993) 

Simultaneous CN and 
gas phase measurements 

Open ocean and coastal 
surface; coastal 
boundary layer i 

Quinn et al (Quinn et al., 
1993) 


Marine boundary layer 

Covert et al. (Covert et 
al., 1992) 


Lower continental 
troposphere 

Radke and Hobbs 
(Radke and Hobbs, 
1991) 

Near small cumulus 
clouds 

Marine free troposphere 

Hegg et al. (Hegg et al., 
1990; Hegg et al., 1991) 

In clear air and near 
stratiform clouds 

Marine boundary layer 

Hoppel et al (Hoppel et 
al., 1990) 

Tropical Atlantic ocean 

Marine boundary layer 

Hoppel and Frick 
(Hoppel and Frick, 
1990) 

Tropical Pacific ocean 

Background continental 
boundary layer 

Marti (Marti, 1990) 

Mountain station, above 
surface inversion layer 
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Abstract 

Sulfates make up a major part of the submicron atmospheric aerosol across a wide range of 
environments. The sulfate aerosol often consists of aqueous sulfuric acid solutions or acidic 
solutions neutralized with sodium, calcium or ammonium. Due to the very low pressures involved 
(below 10“* Pa, or 10*® torr), few measurements of H 2 SO 4 vapor pressure have been made for 

sulfuric acid in the temperature and concentration ranges of atmospheric interest, and no such 
measurements appear to have been made for sulfuric acid solutions neutralized with ammonia. 

This work presents measurements of H 2 S0 4 vapor pressure for aqueous sulfuric acid solutions 
between 55 and 77 wt% H 2 S0 4 (corresponding to about 5 to 25% relative humidity), pure 
ammonium sulfate solids at low humidities, and partially neutralized solutions with [NH 4 MSO^ 
ratios between 0.13 and 1.0. The vapor pressure data collected over sulfuric acid soluuons 
generally agree with the predictions of Ayers, et al [1980], although positive deviation was 
observed for the more dilute solutions. The good agreement between this and previous 
measurements by absolute techniques suggests that the evaporative coefficient for the H 2 S0 4 -H 2 0 
system is near unity. H 2 S0 4 vapor pressures over solid ammonium sulfate were measured 

between 27 and 60C; the data were fit to log p = A/T + B, with A = -2671 ± 369 and B = -1.287 ± 
0.093. The H 2 S0 4 vapor pressures of mixed H 2 S 0 4 -H 2 O-(NH 4 ) 2 SO 4 solutions dropped 
significantly as the [NH^tSO/l ratio exceeded 0.5. The results suggest that ammonia could 

very effectively stabilize molecular clusters of sulfuric acid and water in the atmosphere against 
evaporation, leading to rates of new particle formation higher than predicted by binary H 2 S0 4 -H 2 0 

theory. 



Introduction 


The background atmospheric aerosol varies in composition geographically and temporally. 
Aerosols are commonly found to contain sulfates, including sulfuric acid (H 2 SO 4 ), ammonium 

sulfate ((NH^jSO^ and partially neutralized species such as ammonium bisulfate (NH 4 HSO 4 ). 
Stratospheric aerosols are usually sulfuric acid solutions in the 75 wt% H 2 SO 4 range [Junge, 
1954, 1963], while the tropospheric sulfate aerosol may occur as sulfuric acid or as neutralized 
solutions [Tworaey, 1971; Meszaros and Vissy, 1974; Weiss et al., 1977]. Conventional theory 
holds that these particles nucleate from sulfuric acid and water vapor [Raes and Van Dingenen, 
1992; Easter and Peters, 1994; Raes, 1995] and then undergo neutralization by gas phase 
ammonia. The resulting aerosol consists of a solution of sulfuric acid, water, and ammonium 
sulfate in varying proportions depending on growth history. 

Because of the important roles that aerosols play in the Earth's radiative balance and atmospheric 
chemistry, the nucleation and growth of the sulfate aerosol must be included in global atmospheric 
models. Key to the success of these aerosol models is a thorough understanding of the 
thermodynamics of the ^SC^JkO-CNH^SC^ system. In particular, the values used for 

equilibrium vapor pressures of the system components will profoundly affect predicted aerosol 
nucleation rates. Hence it is important to establish vapor pressure relations for all the components 
of this system within the temperature and solution strength range likely to be found in the 
troposphere. 

Measurements of water vapor pressure over sulfuric acid solutions have been available for decades 
[Grollman and Frazier, 1925; Homung and Giauque, 1955; Giauque et al., I960], in addition to 
values calculated from thermodynamic quantities [Gmitro and Vermeulen, 1964]. The interested 
reader is referred to the comprehensive review of water activity measurements by Bolsaitis and 



Elliot [1990]. Water activities and vapor pressures have recently become available for the H 2 S0 4 - 
H 2 0-(NH 4 ) 2 S0 4 system [Kim et al., 1994; Spann and Richardson, 1985] and for pure aqueous 
(NH 4 ) 2 S0 4 [Apelblat, 1993]. These studies all addressed water, the system component that is by 

far the majority of vapor phase. The vapor pressure of ammonia in these systems is lower by 
many orders of magnitude and is more difficult to measure. Scott and Cattell [1979] succeeded m 
doing so for pure ammonium sulfate between 45 and 180° C; more recently, ammonia vapor 
pressures over ammonium sulfate/acid solutions were reported by Koutrakis et al [1993]. 

The vapor pressure of H 2 S0 4 also lies far lower than that of water for the H 2 S0 4 -H 2 O- 
(NH 4 ) 2 S0 4 system, and is concomitantly difficult to measure. Data on H 2 S0 4 vapor pressures 
over mixed H 2 S0 4 -H 2 0-(NH 4 ) 2 S0 4 solutions appear to be entirely missing from the literature. 
Vapor pressure measurements over unneutralized, concentrated H 2 S0 4 -H 2 0 solutions have been 
accomplished through several ingenious techniques, including radioisotope labeling [Roedel, 

1979], cold trapping of acid vapor [Ayers et al., 1980], and particle shrinkage during levitation 
[Richardson et al., 1986]. These measurements, however, are generally limited to solutions much 
more concentrated (greater than 90% acid by weight) or temperatures much higher (over 50°C) than 
those normally encountered in the atmosphere. In general, atmospheric scientists have made use 
not of the experimental data but of thermodynamic derivations of the H 2 S0 4 vapor pressure over 

sulfuric acid. Older literature cites the derivation of Gmitro and Vermeulen [1964], who calculated 
H 2 S0 4 vapor pressures of sulfuric acid-water solutions starting from pure component 

thermodynamic values. Later work [Verhoff and Banchero, 1972; Banchero and Verhoff, 1975] 
noted that the Gmitro and Vermeulen formulation was inconsistent with previous experimental 
data, and suggested this may be due to uncertainties in the thermodynamic quantities on which the 
formulation was based. Ayers et al [1980] fit their experimental data to the same form as Gmitro 

and Vermeulen, resulting in a formulation that was consistent with the Roedel [1979] 
measurement, but predicted H 2 S0 4 vapor pressures lower than the Gmitro and Vermeulen values 

by over an order of magnitude. The Ayers et al equation for H 2 S0 4 vapor pressures over sulfuric 



acid solutions has been used extensively in sulfate aerosol nucleation modeling, although it is 
based on data obtained at high temperatures (65 to 172 C) and nearly pure acid solutions (98.01 
wt %). To make useful predictions for the troposphere requires an extrapolation in pressure over 

several orders of magnitude. Experimental verification of the Ayers et al fit at lower temperatures 
and solution strengths would be desirable. In addition, there is a clear need for data on H 2 S0 4 

vapor pressures over mixed H 2 S0 4 -H 2 0-(NH 4 ) 2 S0 4 solutions, which probably make up much 
of the submicron atmospheric aerosol. 

This paper presents measurements of H 2 S0 4 vapor pressures over aqueous sulfuric acid 
solutions, solid ammonium sulfate, and mixed H 2 S0 4 -H 2 0-(NH 4 ) 2 S0 4 solutions. Vapor 

pressure measurements were made in temperature and solution composition ranges close to those 

found in the atmosphere, and should be relevant to two current questions regarding the 
atmospheric aerosol. Firstly, which of the H 2 S0 4 vapor pressure formulations found in the 

literature should be used in nucleation models? Secondly, how does H 2 S0 4 vapor pressure 

change as the solution is progressively neutralized with the ammonia? The addition of ammonia 
would presumably cause H 2 S0 4 vapor pressure to drop, thereby stabilizing the solution against 

evaporation. Hence the quantitative behavior of the ternary solution must be known to properly 
model the formation and growth of atmospheric sulfate aerosols. 

Apparatus and Procedure 

Two different experimental approaches were used to measure sulfuric acid vapor pressures. Both 

techniques involved measuring the rates of evaporation of sub micron solid particles and liquid 
droplets into a carrier gas that was initially free of H 2 S0 4 . Vapor pressure data on H 2 S0 4 -H 2 0- 
(NH 4 ) 2 S0 4 solutions, and most of the data for sulfuric acid-water solutions, were obtained by 
measuring the rate at which gas phase H 2 S0 4 evaporating from 0.25 pm particles accumulated in a 
short residence time (2-6 second) flow system. Gas phase measurements were made with a 
chemical ionization mass spectrometer (CIMS), which is described in detail elsewhere [Eisele, 



1986, 1988, 1993; Eisele and Tanner, 1991]. The CIMS system enabled detection of H 2 S0 4 

vapor down to below 10 5 molecules per cnA These measurements of initial evaporation rate 
yielded equilibrium H 2 S0 4 vapor pressure through application of mass transfer theory, as will be 

discussed below. 

Vapor pressures for sulfuric acid-water solutions were also obtained with a tandem differential 
mobility analyzer (TDMA) [Rader and McMurry, 1986]. This technique measured the initial and 
final sizes of monodisperse sulfuric acid solution droplets which flowed through a conditioning 
tube for 30 to 250 seconds. Mass transfer theory was used to calculate equilibrium vapor pressure 
from particle shrinkage. The TDMA technique has been used to measure vapor pressures of 
organic compounds [Rader and McMurry, 1987; Tao and McMurry, 1989]; the approach was first 
applied to the H 2 S0 4 -H 2 0 system by Richert [1991]. 

Both the CIMS and TDMA experiments applied the appropriate mass transfer relations. Since 
mass transfer between particle and gas phase is a function of particle radius (see below), 
monodispersity of sample particles is required. The aerosols (solid crystalline particles or liquid 
droplets) used in the both sets of experiments were size selected by a differential mobility analyzer 
(DMA) [Liu and Pui, 1974]. While the TDMA measurements required significant particle 
shrinkage (in excess of 10%), the mass spectrometric technique required only enough evaporation 
to yield detectable gas phase concentrations; the final sizes of the particles in the CIMS experiments 
were virtually unchanged. 

CIMS measurements. The flow system and detection apparatus are shown in Figure 1. Major 
components included gas phase detection system (the CIMS), an atmospheric pressure flow tube 
system with movable injector, and a particle generation and size selection system, which provided 
monodisperse aerosols of variable composition. 



Liquid droplets of aqueous sulfuric acid or sulfuric acid-ammonium sulfate mixtures were formed 
by atomizing bulk solutions with the desired [NH 4 + ]:[S0 4 = ] molar ratios. The droplets were 

passed through a diffusion dryer containing silica gel, resulting in a somewhat more concentrated 
solution droplet The particles were then imparted with an electrostatic charge by a polonium 
source and size selected with a DMA. A condensation nucleus counter (CNC) (Model 3760, TSI, 
St Paul, MN) was used to monitor output particle concentration. Filtered dry nitrogen was used 
as the carrier gas for aerosol formation and maintained at a flow rate of 1.0 liters per minute. 

The flow system consisted of a temperature controlled Pyrex tube, 80 cm in length by 8 cm in 
diameter, and a movable 2.5 cm diameter glass injector. Monodisperse aerosol from the DMA was 
mixed with particle free nitrogen gas within the injector, which was designed to ensure turbulent 
mixing of aerosol, and flow gas through use of a constriction that raised the flow Reynolds number 
to around 10 4 . A wire mesh covered the flared opening of the injector to distribute the well mixed 
aerosol into the full bore of the flow tube. 

Following particle-gas mixing in the injector, the liquid water content of the droplets very rapidly 
equilibrated with the water partial pressure in the nitrogen carrier gas [Richert, 1991], subject to 
any Kelvin curvature correction. Hence for liquid particles, the flow tube relative humidity 
controlled droplet water content and, by extension, the H 2 S 0 4 concentration of the droplet 

Droplet water vapor pressures were calculated by reference to thermodynamic tables of water vapor 
pressures for aqueous sulfuric acid solutions [Gmitro and Vermeulen, 1964]. The relative 
humidity of the nitrogen flow was controlled by passing some of the gas over a distilled water 
reservoir. Flow gas RH was varied from near zero to about 25% in this manner, and was 
measured by a dew point hygrometer (General Eastern Instruments Model Hygro M-3) at the flow 
tube exit. 


Sulfuric acid-water solution droplets were generated in the 55 to 77 wt% (18 to 36 mol%) range, 



corresponding to the relative humidity range quoted above. Compositions of the mixed H 2 S0 4 - 
H 2 0-(NH 4 ) 2 S0 4 solutions, expressed as the molar [NH 4 + ]:[S0 4 ] ratio, ranged from zero (i.e., 
unneutralized sulfuric acid-water solutions) to 2.0 (pure ammonium sulfate). These measurements 
were performed at relative humidities below 15%, which is weU below both the deliquescence 
point (80% RH) and the crystallization point (39%) of ammonium sulfate [Tang and Munkelwitz, 
1984] and the deliquescence point (39% RH) of ammonium bisulfate ([NH 4 ]:[S0 4 = ] ratio - 1) 
[Tang and Munkelwitz, 1977]. Hence mixed H 2 S 0 4 -H 2 0 -(NH 4 ) 2 S 0 4 particles with 
composition ratio 1.0 to 2.0 were almost certainly solid, crystalline particles; all other aerosols 
measured in the experiment were liquid solution droplets. 


Flow tube temperature was maintained with a circulating water bath; the carrier gas stream was 
passed through a heat exchanger at bath temperature before passing into the injector. Gas 
temperature was checked at the flow tube inlet and exit with calibrated thermocouples. Particle 
residence time in the flow tube was determined by turning the DMA off and on at various injector 
positions, resulting in a rapid decay or rise in particle concentration at the CNC located at the flow 
tube exit By varying the injector position and establishing the delay between DMA transition and 
measured changes in particle counts, flow tube residence times for each injector position were 
established. These times, ranging between 2 and 6 seconds, were equal (within measurement 
error) to those predicted by laminar flow. 

The flow tube was directly coupled to the CIMS system. The mass spectrometer measured count 
rates at mass 97 (the HS0 4 ‘ ion) and mass 62 (the NO 3 - ion). The (97)/(62) count ratio, 
multiplied by a predetermined calibration factor, yielded the total H 2 S0 4 concentration in the gas 
flow stream with an absolute accuracy of about 30 to 35%. Instrument background was obtained 
by measuring the count rates from several masses, such as mass 10 and 106, where no signal was 
expected to occur. Calibration of the mass spectrometer is described by [Tanner and Eisele, 1995] 
A CNC (TSI 3760) measured aerosol concentration at the flow tube exit. Measurements of both 



particle and gas phase concentration were typically integrated over 30 second intervals. 

TDMA measurements. The TDMA technique has been discussed in detail by [Rader and 
McMurry, 1986], [Tao and McMurry, 1989], and [Richert, 1991]. The TDMA apparatus is 
shown in Figure 2. Cleaned, filtered air flow from an air purifier (Aadco Model 737) was split and 
bubbled separately through reservoirs of distilled water and fuming sulfuric acid, the latter 
becoming enriched in sulfur trioxide, SO 3 . The flows were then turbulently combined in a small 
mixing volume, forming droplets of concentrated sulfuric acid- water solution. These particles 
were electrostatically charged and passed through a DMA to produce a monodisperse aerosol; sizes 
between 0.05 and 0.3 |im diameter were used in the experiment The ability of the DMA to 
accurately size particles above about 0. 1 |im is limited by multiple charging. For this portion of the 
experiments, particles were charged with the low ion concentration, bipolar device described by 
Gupta and McMurry [ 1 989] to reduce the effects of multiple charging. 

•4 

The 1.5 1pm aerosol flow was then mixed with 8.5 1pm of purified air and passed into a 3 m 
conditioning tube, 10 cm in diameter. The conditioning tube temperature was controlled by a 
circulating bath. The 10 1pm tube flow was laminar, however, residence time measurements, 
similar to those described above, yielded particle residence times somewhat longer than predicted 
by laminar flow. As with the flow tube experiments described above, the conditioning tube flow 
gas was humidified as desired, and droplet composition was determined assuming the water vapor 
pressure of the droplets to be in equilibrium with the flow gas. 

Theory 

Transition regime mass transfer. Equilibrium vapor pressure of H 2 SO 4 was determined from the 
experimentally measured quantities by mass transfer theory. With particle radius set at 0.125 [im, 
the Knudsen number Kn of the system (= molecular mean free path/particle radius) is about 0.25, 
placing mass transfer in the transition regime. Several transition regime mass transfer expressions 



10 



where D vg is the diffusivity of the vapor species, M the vapor molecular weight, D p the particle 
diameter, p vs the equilibrium vapor pressure and p^ the partial pressure of the vapor far from the 
particle. Two expressions for F(Kn) were considered in this work, those of Fuchs and Sutugin 
[1970] 


F(Kn) = — = 1 + Kn 


1.333+.71Kn 
l + Kn~ J 


and Bademosi and Liu [1971 a,b] 

(4) 

The ratio Jc/Jk is that of continuum to free molecular mass transfer expressions, which can be 
written as [Davis and Ray, 1978] 

J c _ 4j3Kn _ 4D vg 
J k E Eac 



( 5 ) 



where |5 is a constant obtainable from Enskog-Chapman theory for molecular diffusivity, a is 
particle radius and c is the mean molecular velocity. In expressions (3) and (5), E is the 
evaporative coefficient, which will be addressed in detail below. Substituting (5) in (4) yields a 
function of E and Kn analogous to equation (3): 



The transition regime corrections given by equations (3) and (6) agreed to within about 10%. 

While the Fuchs-Sutugin expression is widely used for transition mass transfer, Davis and Ray 
[1978] found that the Bademosi-Liu relation better fit their experimental data for Kn > 0. 1. For 
this reason equation (6) was employed in this work. 

CIMS Measurements. Each experimental run with the CIMS system measured the accumulation of 
gas phase HjSC^ molecules during evaporation from an ensemble of 1000-3000 particles cm' * 

The total mass flux from n p monodisperse particles of diameter Dp is 

^ = n p nD p 2 J = n p 7tD p 2 F(Kn)J c (7) 

This mass flux evolves a (very low) partial pressure given adequately by the ideal gas law 

dm _ MV dp .q, 

dt~ kT dt K } 


Combining equations 4, 7, and 8, we get 



^ = 2 F( Kn)NnD p D vg ( p VJ - pjt )) (9) 


where N = rij/V, the number density of aerosols. The experimentally measured quantities in (9) 
are aerosol concentration N (fixed), H 2 SO 4 partial pressure P^fat ^t^^ ), and Dp (fixed). F(Kn) 
is a constant for fixed D p . Equation 9 may be integrated to yield an expression for equilibrium 
vapor pressure [Camraenga, 1980] 

p _ Pm ~ Po exp(-Cr) ( 10 ) 

” 1 - exp(-Cr) 

where C is the collection of constants in (9) 

C = 2F( Kn)NnD p D vgt 

p m is the H 2 S0 4 pressure measured at the end of the flow tube, and p 0 is the very small initial 
concentration of H 2 S 0 4 brought into the injector along with the aerosol flow (a result of the sulfate 

aerosol beginning to evaporate before reaching the flow tube entrance). This initial pressure was 
estimated to be less than 10 ' 9 Pa, based on expected evaporation rates of the concentrated acid 
solution droplets before they were mixed with humidified flow gas. The ratio p</p m was generally 

below 10%, keeping any error introduced by this estimate tolerably smalL In general, equation 10 
should include a correction for the elevation of droplet vapor pressures due to the Kelvin effect 
For the droplets used in this portion of the study (D p = 0.25 pin) the Kelvin correction amounted 

to less than a tenth of a percent and was ignored. 

TDMA Measurements. After evaporation times of 30 to 250 seconds, the sulfuric acid solution 
droplets were sampled from the TDMA conditioning tube and their final size determined by a 
second DMA The change in particle diameter D p is related to mass flux density J by 
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d£r _ l 

dt p 


( 11 ) 


where p is solution droplet density. Using the Bademosi-Liu transition regime mass transfer 
equation (3) and incorporating a Kelvin curvature term (droplets used in TDMA experiments were 
small enough to require the Kelvin correction), equation 1 1 can be rewritten and integrated [Tao 
and McMurry, 1989]: 


P v,0 - P- + 


pRT 

4D vt Mf 


f O.J, 

r -4eM y 

J*,7 exp 

[pRTDj 


D dD 

P P 


( 12 ) 


where the integration limits are initial to final droplet diameter. Equation 9 was solved numerically, 
using measurements of Dj, D f , residence time t and tube temperature, and estimates of Poo, the 

"back pressure" of H2SO4 that accumulated in the conditioning tube over time. Since the solution 
droplet composition was fixed by carrier gas relative humidity, any H2S0 4 evaporation would be 

accompanied by associated water to keep the droplet composition constant To accurately account 
for particle shrinkage, both the H 2 S0 4 and H 2 0 mass fluxes were calculated separately. 

The mass transfer expressions (3) and (6) contain the evaporation coefficient E. For simple 
liquids, E is taken to be equal to the coefficient of condensation otc. Either coefficient must be 
determined experimentally for individual compounds; E (otc) is defined as the measured rate of 
evaporation (condensation) divided by the rate theoretically predicted with E (otc) set equal to 
unity. While the value of otc for sulfuric acid solutions has been reported as much less than unity, 
between 0.02 and 0.09 [Van Dingenen and Raes, 1991], there is a considerable body of work that 
suggests that for single component liquids, E is unity or very close to unity [Mozurkewich, 1986]; 
Cammenga, 1980; Pound, 1972; Maa, 1967] and that experimental results to the contrary may be 
due to sample cooling and contamination. The systems described in this work are multicomponent 



solutions, so these arguments may have limited applicability. However, Mozurkewich [1986] has 
argued that otc for any polar species on an aqueous surface should be near unity. This would 
imply that E for a sulfuric acid solution or an ammonium sulfate-sulfuric acid mixture should be 
near unity, even though they are not pure liquids. Due to the undetermined state of E for the 
sulfate systems examined in these experiments, the vapor pressure data presented in this work have 
been calculated with E assumed to be unity. The validity of such an assumption may be judged by 
comparing the results of this work with extrapolations of measurements at higher temperatures and 
solution strengths [Ayers et al., 1980; Roedel, 1979] which did not rely on mass transfer 
considerations. 

Results 

The vapor pressures of HjSC^ over aqueous sulfuric acid solutions at 25, 30 and 35°C are shown 
in Figures 3a-c. The data in each figure are for a single temperature, with solution concentration 
varied from about 55 to 77 wt%. Figure 3a shows vapor pressures measured with the CIMS only. 
Figures 3b and 3c include both CIMS and TDMA derived data. In all cases, the data fall closer to 
the predictions of Ayers et al than to the earlier Gmitro and Vermeulen work. However, the data 
agree with the Ayers et al formulation only for solution concentrations above 70 wt%. Below this 
solution strength the vapor pressure data appear to be up to 5 times higher than the predictions of 
Ayers et al, although still well below those of Gmitro and Vermeulen. Reasons for the discrepancy 
are explored below in Discussion. 

Vapor pressure measurements obtained with the TDMA generally support those from the CIMS 
experiments (Figures 3b, 3c). However, the figures show the high scatter that the TDMA 
technique was prone to, scatter well in excess of random error. Both the present experiment and 
previous work [Richert, 1991] found that calculated H 2 S0 4 vapor pressures often depended on 

particle residence time within the TDMA. Vapor pressures derived after shorter evaporation times 



(t< 30 seconds) were usually nearer to those predicted by Gmitro and Vermeulen but dropped as 
evaporation time increased over 60 seconds, reaching near zero after the full 250 second tube 
residence time. TDMA results such as these suggested the presence of a droplet contaminant that 
lowered the l^SO^ vapor pressure, such as ammonia gas, the effects of which would become 

more pronounced with dropping particle size. The TDMA system and associated plumbing was 
mechanically and thermally cleaned and examined for possible contamination problems. These 
efforts reduced data scatter somewhat and yielded the data shown in Figures 3b and 3c. However, 
the performance of the TDMA with the HjSO^-^O system did not approach that of the CIMS, as 

is evident in the Figures. The better repeatability and accuracy of the CIMS system led to reliance 
on this technique for the bulk of the ^SO^-water measurements and all the measurements of 

mixed systems. 

Measured H 2 SO 4 vapor pressures of pure (NH^SC^ are shown in Figure 4 for the temperature 

range 28 to 60° C The relative humidity for all measurements was kept below about 6 %. The data 
are plotted as log pressure (in Pa) versus inverse temperature (K*^). Although possessing 
substantial uncertainties, the points form a convincing regression line of the Clausius-Clayperon 
form 

log p = ~ + B (13) 

with A = -2670 ± 369 and B = -1.29 ± 0.0925. The extremely low vapor pressures measured in 
this temperature range were near the lower detection limit of the apparatus, as evidenced by the 
large error bars on the data. 

As noted above, the (NH^SO^ particles were most likely solid, since the relative humidity of the 
system was well under the crystallization point of ammonium sulfate. The particles’ vapor 
pressure-temperature relationship is unlikely to change with relative humidity as long as they 



remain in the solid phase. Starting as such, ammonium sulfate aerosols can be expected to stay 
solid until the ambient relative humidity reaches the deliquescence point Hence the relationship in 

equation 13 should hold for atmospheric humidities below the deliquescence point The only 
prediction of H2SO4 vapor pressure over ammonium sulfate found in the literature is that of Scott 

and Cattell [1979], who calculate a vapor pressure of about 10'13 Pa for nearly pure (NH^SC^ 
at 25 °C. This is two orders of magnitude lower than given by the experimental fit equation above. 

In Figure 5 are shown vapor pressure data from mixed solutions of sulfuric acid, ammonium 
sulfate and water, collected across a range of relative humidity while keeping the temperature fixed 
at 30° C. The data are grouped by [NH^tlSO^ ratio, hereafter called the ionic ratio. The vapor 

pressures measured over solutions with low ionic ratios (below 0.2) were only slightly below 
those measured for unneutralized sulfuric acid solutions, especially at higher humidities. Solutions 
with ionic ratio = 0.33 and lower also possessed vapor pressures close to the sulfuric acid values. 
It was not until the ionic ratio increased over 0.5 that H2SO4 vapor pressure was appreciably 
depressed vis-a-vis the unneutralized acid. Particles with ratio = 0.8 had H2SO4 vapor pressures 

less than one quarter of those for sulfuric acid solutions at 15% relative humidity; at 5% relative 
humidity they were over an order of magnitude lower than their unneutralized acid counterparts. 
Particles with ionic ratio of 1.0 (probably ammonium bisulfate, NH 4 HS0 4 ) exhibited H 2 S0 4 

vapor pressures in the mid 10‘ 10 Pa range that were nearly constant between 5 and 15% relative 
humidity. Data take n at this temperature for ammonium sulfate, although sparse, also show no 
apparent change with relative humidity, as would be expected for the solid below its deliquescence 
point. The vapor pressures for the pure ammonium sulfate solid appears to be below 10* 10 Pa at 
these relative humidities. 

Discussion 

The H 2 S0 4 vapor pressure measurements presented in this paper depend in part on two key 
parameters, the evaporation coefficient E and the H2S0 4 diffusion coefficient D. The values 



selected for these parameters affect the outcome of equilibrium vapor pressure calculations. As 
addressed above in the Theory section, this work assumed an evaporative coefficient of unity. The 
vapor pressures derived with this assumption, shown in figures 3a-c, are in substantial agreement 
with the extrapolation of the Ayers et al fit line in the higher solution concentration range. The 
latter is a fit to vapor pressure measurements obtained at higher temperatures and solution 
strengths; these measurements were absolute, i.e., they did not rely on mass transfer assumptions. 

Thus if the Ayers et al results can be extrapolated, the present work offers evidence that E is close 
to unity for the H2SO4-H2O system. There is a clear need for further experiments in this area. 

Vapor pressures presented in this paper were derived using value of the diffusion coefficient D for 
H 2 S0 4 in nitrogen/air calculated from Enskog-Chapman theory [Present, 1958]. One obstacle to 

using this approach is the selection of a suitable diameter for the collision volume of the highly 
aspherical H 2 S0 4 molecule; small variations in the choice result in large changes to the calculated 
D vg . In addition, conditions in this experiment were such that the gas phase H 2 S0 4 molecules 

were hydrated with one or more water molecules [Heist and Reiss, 1974]. For this work, a 
collision diameter for H 2 S0 4 was chosen as that of a hydrated H 2 S0 4 cluster from the work of 
Heist and Reiss [1974]. The value of D vg so obtained was 0.103 ± 0.01 cm 2 sec' 1 at 25C, 
ranging up to 0.1 10 at 60C. This compares with D suggested (but not explicitly derived) by 
Roedel [1979], 0.08 ± 0.02 cm^ sec -1 at 23C, which has become widely used in the atmospheric 
science literature. 

The vapor pressures of sulfuric acid solutions (Figures 3a-c) are largely consistent with the 
predictions of Ayers et al [1980] at solution concentrations greater than 70 wt%. In this range the 
data clearly support the 1980 Ayers et al formulation for H 2 S0 4 vapor pressure over the Gmitro 

and Vermeulen predictions. Below 70 wt% however the data positively diverge from the Ayers et 
al line, moving closer to, but still substantially below, its older counterpart. Two possible causes 
for such an effect are droplet cooling, which should lower the vapor pressure of the more volatile 



(higher wt%) droplets, and particle shrinkage, which would give erroneous calculated vapor 
pressures through its effect on particle surface area and thus mass flux density. Final droplet 
temperatures were estimated by means of a coupled mass-heat transfer calculation [Vesala, 1991]. 
Thermal conductivity and specific enthalpy were assumed constant over the relatively narrow 
temperature range examined in this study. The results indicated a negligible cooling effect even for 
the most concentrated (77 wt%) solution droplets. The very low decrease in droplet temperature 
(less than 10"^ K) was due to the minute amounts of liquid phase material that volatilized in the 
short tube residence time. These short evaporation times likewise kept calculated particle shrinkage 
to be less than 0. 1 % of initial volume. It seems then that the divergence of H 2 S0 4 vapor pressures 

with more dilute solutions, shown in Figure 3, may be real. 

Recent field measurements have shown a low but nonzero nocturnal partial pressure of H2SO4 
[Eisele,1993]. This is an unexpected result since the production of gas phase H2SO4 relies on the 
oxidation of SO2 by OH; the latter, being in rapid photochemical equilibrium, vanishes quickly 
after sunset, thereby turning off the major H2SO4 source. In the absence of a photochemical 
source, the H2SO4 may be present due to the vapor pressure of the background aerosols. Eisele 
and co-workers [1993] measured roughly 4x10'^ Pa H2SO4 at a clean marine site in Washington 

state and about 4x10'^ p 4 at a station in the mountains of Colorado. Although both sites 

probably had relative humidities in excess of those used in this present experiment, the data in 
Figure 4 may be extrapolated for use as a rough guide to H2SO4 vapor pressures for moister 

conditions. The H2SO4 partial pressures measured in at these field sites may be more consistent 

with a partially neutralized sulfate droplet; an extrapolation of the data in Figure 5 to higher relative 
humidities suggest an aerosol with [NH^"*"] : [ S 0^ = ] ratio of 0.8 or less. 

Data for H 2 S0 4 -H 2 0-(NH4) 2 S0 4 mixed solutions (Figure 5) are also subject to the evaporation 

coefficient concerns outlined above. Regardless of their absolute values, however, the vapor 
pressures over H 2 S0 4 -H 2 0-(NH 4 ) 2 S0 4 mixed solutions convey important information from the 



relative placement of the curves. As Figure 5 shows, addition of the ammonium ion to a sulfuric 
acid-water solution begins to significantly depress H 2 S0 4 vapor pressure as the ammonium to 

sulfate ion ratio exceeds about 0.5. These data were obtained using sub micron droplets, a 
collection of about 10 7 molecules. If clusters of 1 to 10 H 2 S0 4 and H2O molecules behave in an 

analogous manner. Figure 5 suggests that a few ammonium ions could substantially stabilize such 
molecular clusters. The result would be smaller critical cluster sizes and higher observed particle 
formation and growth rates in the presence of ammonia. Since levels of atmospheric ammonia 
have been found to be at least an order of magnitude greater than measured H 2 S0 4 concentrations' 

at continental sites [Tanner and Eisele, 1991; Langford and Fehsenfeld, 1992], it is reasonable to 
suppose that ammonia may play a significant role in new particle formation at these sites. Recent 
field measurements by the authors [Weber et aL, 1995b] indicate that new particle formation rates 
frequently exceeded those predicted by binary nucleation theory at a remote continental site. These 
measurements offer strong evidence that some other substance may be involved in stabilizing 

«4 

H 2 S0 4 -H 2 0 clusters. The present work provides quantitative support for atmospheric ammonia as 
one such stabilizing agent 
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Figure Captions. 


Figure 1. Chemical ionization mass spectrometer and flow tube apparatus used in vapor pressure 
measurements. Aerosol flows are given with solid arrows; sheath or carrier gas flows are depicted 
with open arrow symbols. 

Figure 2. Tandem differential mobility analyzer apparatus used in vapor pressure measurements. 

Figure 3. H 2 S0 4 vapor pressures (in Pascals) measured over sulfuric acid solutions of varying 
concentration. Labeled lines show two vapor pressure formulations from the literature. 

3a. H 2 S0 4 vapor pressures at 25 °C. 

3b. H 2 S0 4 vapor pressures at 30 °C. Both CIMS and TDMA derived data 
are included. 

3c. H 2 S0 4 vapor pressures at 35 °C. Both CIMS and TDMA derived data 
are included. 

Figure 4. H 2 S0 4 vapor pressures over samples of ammonium sulfate at low (<5%) relative 
humidity. 

Figure 5. H 2 S0 4 vapor pressures of partially to fully neutralized sulfuric acid. Data are grouped 
by the ionic ratio (see text). 
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